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Abstract. We introduce the concept of "D-operators associated to a sequence 
of polynomials and an algebra A of operators acting in the linear space 

of polynomials. In this paper, we show that this concept is a powerful tool 
to generate families of orthogonal polynomials which are eigenfunctions of a 
higher order difference or differential operator. Indeed, given a classical dis- 
crete family (pn)n of orthogonal polynomials (Charlier, Meixner, Krawtchouk 
or Hahn), we form a new sequence of polynomials {qn)n by considering a linear 
combination of two consecutive pn- Qn = Pn + finPn—i, fin G K- Using the 
concept of X>-operator, we determine the structure of the sequence {f)n)n in 
order that the polynomials {q„)„ are common eigenfunctions of a higher order 
difference operator. In addition, we generate sequences {l3n)n for which the 
polynomials {q„)„ are also orthogonal with respect to a measure. The same 
approach is applied to the classical families of Laguerrc and Jacobi polynomi- 
als. 
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1. Introduction and results 

The issue of orthogonal polynomials {pn)n which are common eigenfunctions 
of a difference operator goes back at least for one century and a half. The first 
example was introduced by Chebyshev in 1858. Along the first decades of the XX 
century, some other examples appeared and they are associated with the names of 
Charlier, Meixner, Krawtchouk and Hahn. In 1941, O.T. Lancaster [31] classified all 
orthogonal polynomials on the real line satisfying second order difference equations 
of the form 

(tAVp„ -f rVpn = A„p„, n > 0, 
where a and r arc polynomials of degree at most 2 and 1, respectively (independent 
of n), and A and V denote the first order difference operators: 

(1.1) A(/) = /(x + 1) - /(a:), V{f)^ f{x)- f{x-l) 

(that is, the polynomials {pn)n are common eigenfunctions of the second order dif- 
ference operator crAV -I- rV). If we consider orthogonal polynomials with respect 
to positive measures (on the real line), the only solutions are the polynomials of 
Charlier, Meixner, Krawtchouk and Hahn. For similarity with the families of Her- 
mite, Laguerre and Jacobi, they are usually called classical discrete polynomials. 
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The adjective discrete means that each of these famihes is orthogonal with respect 
to a discrete measure. More precisely, the discrete orthogonalizing weights for these 
families of polynomials are {N e N \ {0}): 



(1.2) Charlier y^—r^x, < a; 

(1.3) Meixner ^ Q^r(x + c) ^^^ < a < 1, < c: 

N-l ^ 

(1.4) Krawtchouk t^tttt n^^^ < a; 

■^-^ 1 (A* — x)x\ 

a — x)T{: 
T{N -x)xl 



x=0 

N-l, 



x=0 



However, almost nothing was known about orthogonal polynomials which are 
common eigenfunctions of a higher order difference operator. We will expand a 
finite difference operator using the basis formed by the shift operators 5;, Z e Z, 
defined by 

"difih] (1.6) 5i{f)^fix + l) 

(this basis is equivalent to the basis formed by the powers of the difference operators 
A and V defined by (jl.ip ). Hence we consider finite difference operators of the form 



(1.7) D = J2fi3h s<r,s,reZ. 

l=s 

If fr, fs 7^ 0, the order of D is then r — s. We also say that D has genre (s, r). 

Until very recently no examples of orthogonal polynomials which are common 
eigenfunctions of a higher order difference operator of the form ()1.7|) were known, if 
we except the discrete classical families themselves. What we know for the continuos 
or q cases (differential or g-difference operators, respectively) has seemed to be of 
little help because adding Dirac deltas to the classical discrete families seems not to 
work (see the papers [2^ and [3' by Bavinck, van Haeringen and Koekoek answering, 
in the negative, a question posed by R. Askey in 1991). 

The first examples of orthogonal polynomials {pn)n which are common eigen- 
functions of a difference operator of order bigger than 2 have been introduced by 
the author in [8 . We generate the orthogonalizing measures for these families of 
polynomials multiplying the classical discrete weights by certain variants of the 
annihilator polynomial of a set of numbers. For a finite set of numbers F, this 
annihilator polynomial p_F is defined by 

(1.8) pF{x)^ll{x-j). 

Using this annihilator polynomial (and other of its variants), we have posed in [S] 
a number of conjectures. Two of them are the following. Write np for the number 
of elements of F. 

Conjecture 1. Let p be either the Meixner (II. 3|) or the Hahn weight (|1.5p (in 
this last case we assume N big enough to avoid trivialities). For a finite set F of 
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positive integers, consider the polynomial p-c-F defined by 



where c is the parameter in the Meixner and Hahn weights. Then, the sequence 
(finite or infinite) of orthogonal polynomials with respect to the measure pF,c = 
P-c-fP are common eigenfunctions of a difference operator D of the form 



Conjecture 2. Let p be any of the classical discrete weights of Charlier, Meixner, 
Krawtchouk or Hahn (in this last cases we assume N big enough to avoid triviali- 
ties). For a finite set F of positive integers, assume that the measure pp — PfP has 
associated a sequence {pn)n of orthogonal polynomials (finite or infinite: see the 
preliminaries for more details) , where pF is the polynomial defined by ()1.8|) . Then, 
these orthogonal polynomials {pn)n are common eigenfunctions of a difference op- 
erator D of the form 



The purpose of this paper is to introduce a method to generate sequences of 
orthogonal polynomials which are eigenfunctions of a higher order difference oper- 
ator. The method is based on the concept of X'-operator associated to a sequence 
of polynomials {pn)n and an algebra A of operators acting in the linear space of 
polynomials. Starting from one of the classical discrete families {pn)n of orthog- 
onal polynomials, we form a new sequence of polynomials ((?«)„ by considering a 
linear combination of two consecutive Pn- qn ^ Pn + PnPn-i, Pn G Using the 
concept of ©-operator associated to the family (p„)„ and the algebra of difference 
operators with polynomial coefScients, we determine the structure of the sequence 
(/3n)n in order that the polynomials (g„)„ are commons eigenfunctions of a higher 
order difference operator. We will use this technique to prove the conjectures 1 and 
2 above when F is the finite set formed by the first k consecutive positive integers 
F — {1,2,-- - fc > 1. This method also generates sequences of polynomials 
which are eigenfunctions of higher order differential (see the Appendix of this pa- 
per) or g-difference operators (see [T]). Hence, the concept of 25-operator provides 
an unified approach to construct orthogonal polynomials which are eigenfunctions 
for each one of these three different types of operators (differential, difference or 
g-difference). 

The content of this paper is as follows. 

In Section 3, we consider the key tool of our method: D-operators associated to 
a sequence of polynomials {pn)n (Pn of degree n) and an algebra of operators A 
acting in the linear space of polynomials. We have different types of P-operators. 
The simplest ©-operator is defined from an arbitrary sequence of numbers (e^)™ 
as follows: we consider first the operator C : P — ;> P defined by linearity from 
C{Pn) = £nPn-i- Usiug that dgr(^(p)) < dgr(p) for all polynomial p, we define the 
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operator V by 



(1.9) 



2? : P -> P 

OO 



We then say that P is a P-operator for the sequence (p„)„ and the algebra A if 

For the examples in this paper we consider two algebras A. When the sequence of 
polynomials (p„)„ is one of the classical discrete families (Sections 4, 5, 6 and 7) i.e. 
Charlier, Meixner, Krawtchouk or Hahn polynomials, we consider the algebra A 
formed by all finite order difference operators X]j=s fj^i where fj is a polynomial, 
j = s, • • • , r: 



algdiff (1.10) 



, r, s < r 



When the sequence of polynomials {pn)n is either the Laguerre or Jacobi poly- 
nomials (Appendix), we consider the algebra A formed by all finite order differ- 
ential operators '^j^q fjid/dxy where fj is a polynomial with degree at most j, 
j=0,.-. ,k: 



algdiffl (1.11) 



-4= E/. 



/, eP,dcg(/j) <j,j-0,..,fc,/ceN 



Here it is an example of 2?-operator. Consider the Charlier polynomials (c°)„ 
(see ()4.ip below), a ^ 0, and the algebra A given by ()1.10|) . Then the sequence 
e„ = 1, n > 0, generates a D-operator for the Charlier polynomials and this 
algebra A; indeed, the operator V defined above (II. 9p from the sequence £„ = 1, 
n > 0, and the Charlier polynomials is equal to the Nabla operator (see (jLip ): 
V — \7 ^ A (see Lemma [4. II) . In this paper, we will also display two P-operators 
for Meixner and Krawtchouk polynomials (see Lemmas 15.11 16.1|) , four P-operators 
for Hahn polynomials (see Lemma l7.2p . one for Laguerre polynomials (see Lemma 
18. ip and two for Jacobi polynomials (see Lemma [8. 7p . For P-operators associated 
to g-classical families of orthogonal polynomials see [T] . 

Roughly speaking, our method consists in the following. In addition to the 
polynomials {pn)n and the algebra A, we also have an operator Dp S A for which the 
polynomials (p„)„ are eigenfunctions. We denote (0„)„ the corresponding sequence 
of eigenvalues: Dp{pn) = 6'„p„. To simplify we consider here in the Introduction 
the case when 0„ is linear in n. Assume also that we have identified a sequence 
{£n)n which defines a I?-operator for the sequence of polynomials {pn)n and the 
algebra A. We can then produce a large class of polynomial families {qn)n which 
are eigenfunctions of operators in the algebra A proceeding as follows. Take an 
arbitrary polynomial Q (such that Q{9n) 7^ 0, ?i > 0), and define the sequence of 
numbers (/3„)„>i by 



'"Q(^„-i)' 
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Define next a new sequence of polynomials {qn)n by go = 1 and 



qnint 



(1.12) 



■ PnPn 



n> 1. 



th5. Ich 



becchar 



qnno 



It turns out that the sequence of polynomials {qn)n are eigenfunctions for the oper- 
ator Dg = P{Dp) + 'DQ{Dp), where the polynomial P satisfies P{6„) — P(0„_i) = 
Q{Sn-i) (since {9n)n is lineal in n, P always exists). Since 2? is a I?-operator, 
V ^ A and then Dq £ A as well (see Section 3 for more details). 

Only for a convenient choice of the polynomial Q, the polynomials ((?„)« (I1.12[) are 
also orthogonal with respect to a measure. When the sequence {pn)n is any of the 
classical discrete families of orthogonal polynomials, a very nice symmetry between 
the family {pn)n and the polynomial Q appear. Indeed, if Pn, n > 0, are either the 
Charlier, Meixner or Krawtchouk polynomials, then when we choose the polynomial 
Q in the same family (but with different parameters) the polynomials (g'„)„ turn 
out to be orthogonal as well, while for the Hahn polynomials, the polynomial Q 
has to be chosen in the dual family. For other choices of the polynomial Q, the 
polynomials ((j„)n seem not to be orthogonal with respect to any measure. In a 
subsequent paper, we will show that they enjoy certain orthogonality properties 
which imply that the family (q„)„ has to satisfy certain recurrence relation of order 
at most 2 deg{Q)+3. In the terminology introduced by Duistermaat and Griinbaum 
([B]; see also |T2], [13]), we can say that our method produces bispectral polynomials: 
that is, the family (q'„(a;))„ are common eigenfunctions of two difference operators, 
one in the continuous parameter x and the other in the discrete parameter n. 

Sections 4, 5, 6 and 7 are devoted to show how our method works by applying 
it to the families of classical discrete polynomials. Hence, this paper together with 
the forthcoming one ([T]) show that the concept of 2?-operator provides an unified 
approach to construct orthogonal polynomial which are eigenfunctions of higher 
order differential, difference or g-difference operators, respectively. 

Here it is an example of the kind of results we can get using the concept of 
2?-operator. 

Theorem 1.1. Let Pi be a polynomial of degree k + I, k > 1, and write P2{x) = 
Pi{x~l)~Pi{x) (so that P2 is a polynomial of degree k). We assume that P2{—n) 7^ 
0, n > 0, and define the sequences of numbers 



(1.13) 



In+l 
I3n 



P2{-n), n>0, 
Pi{-n), n>0, 

7n+l 



n > 1, 



In 



and the sequence of polynomials qa = 1, and for n> I 
(1.14) g„(x) =<(x)+ (a;), 

where c°, n > 0, a ^ 0, are the Charlier polynomials (see (^T7p below). Consider 
the second order difference operator Da li4-4^ with respect to which the Charlier 
polynomials are eigenfunctions. Write finally D for the difference operator of order 
2k + 2 and genre (— fc — 1, fc + 1) 

D^Pl{Da)+VP2{Da). 



Then D{qn) = A„q„, n>0. 
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When we choose the polynomial P2 in an appropriate way, the polynomials {qn)n 
()1.14|) are also orthogonal with respect to a measure. As we have pointed out 
above, there is a very nice symmetry in the choice of P2. Indeed, for the Charlier 



case, when we choose ^2(2;) 



, "(x — 1) and assume Cj,°(— n) 7^ 0, n > —1, the 



polynomials {qn)n (|1.14l) are then orthogonal with respect to the measure 

00 x+k+l 



This provides a proof for our Conjecture 2 above for the Charlier polynomials 
and the finite set F = {1, 2, ■ • • , fc}. Indeed, the measure p_FPo, where p_F is the 
polynomial (|1.8|) and pa is the Charlier measure (|1.2|) . is equal to the measure pk,a 
above shifted by A: + 1: pppa = Pk,a{x ~ k ~ \). 

We have added an Appendix devoted to orthogonal polynomials {pn)n which are 
eigenfunctions for a higher order differential operator of the form 



(1.15) 



E 



flPn^ = KPn, n>0, 



where //, Z = 1, • • • , fc, are polynomials independent of n (so that (p„)„ are common 

eigenfunctions of the differential operator D ~ X]/=i // (^)')- This issue goes back 
to the late thirties, when H.L. Krall classified all families of orthogonal polynomials 
which satisfy a fourth order differential equation of the form (jl.isp (see p5t i26) ). 
Besides the classical families of Hermite, Laguerre and Jacobi (satisfying second 
order differential equations) , he found three other families of orthogonal polynomi- 
als satisfying a fourth order differential equation of this type. These polynomials 
are orthogonal with respect to a positive measure which consists of some classical 
weights together with a Dirac delta at the end point(s) of the interval of orthog- 
onality: X[-i,i] + -^(^-i -I- Si) (Legendre type), e""^ -|- MSq (Laguerre type) and 
(1 - a;)"x[o,i] + ^^^0 (Jacobi type). 

Since the eighties a lot of effort has been devoted to this issue. T. H. Koornwinder 
([24]) studied in 1984 orthogonal polynomials with respect to the Jacobi weight 
together with two Dirac deltas at ±1 



jacobid (1-16) 



{l-x)"{l+xf +M6-1+NS1, a,l3>-l. 



As a limit case he found orthogonal polynomials with respect to the Laguerre weight 
together with a Dirac delta at 



MSo, a > -1. 



Some years earlier, L.L. Littlejohn ([SHIM]) had discovered new families satisfying 
sixth order differential equations. They correspond with the case a = /? = in 
(|1.16l) (Krall polynomials; discovered by H.L. Krall but never published) and a — 
1,2 in (jl.l7p (Littlejohn polynomials). A. M. Krall and L. L. Littlejohn did some 
work on the classification of higher order differential equations of the form (|1.15p 
having orthogonal polynomial solutions (for several values of fc > 6) (see p71 153] ). 
J. Koekoek and R. Koekoek showed in 1991 that orthogonal polynomials with 
respect to the weights (|1.17l) satisfy infinite order differential equations, except for 
nonnegative integer values of a for which the order reduce to 2a + A f|20|). Similar 
results were proof by R. Koekoek for the Koornwinder polynomials orthogonal with 
respect to the weight (|1.16p for a = /3 and M = N ([H]), and later on (again with 
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J. Koekoek) for the case M,N > ([21 ). K.H. Kwon and D.W. Lee proved that 
when we consider orthogonal polynomials with respect to a positive and compactly 
supported measure satisfying higher order differential equations, they have to be 
orthogonal with respect to some of the measures (|1.16p ([28), see also [29], [30]). 
A. Zhedanov proposed a technique to construct Krall's polynomials and found a 
{2a + 4)-th order differential equation for the orthogonal polynomials with respect 
to ()1.16|) when a is a nonnegative integer and M — ( 38 ). 

F.A. Griinbaum and L. Haine (et al.) proved that polynomials satisfying fourth 
or higher order differential equations can be generated by applying Darboux trans- 
forms to certain instances of the classical polynomials (see [T2| [14] [TS] [17]). By 
modifying the Darboux process, they also introduced new families of polynomials 
satisfying fourth order differential equations of the form (|1.15|) which also satisfy 
five term recurrence relations, instead of the three term recurrence relation which 
characterizes the orthogonality with respect to a measure (see [13]). P. Iliev refined 
this approach and characterized the algebras of differential operators associated to 
these families of polynomials (see [18], |19j). 

In the Appendix we construct one 2?-operator for Laguerre polynomials and two 
for Jacobi polynomials, and apply our method to recover some of the results com- 
mented in the previous paragraphs. 



2. Preliminaries 

For a linear operator Z) : P — > P and a polynomial P{x) = X)^=o°'i^"'' 
operator P{D) is defined in the usual way P{D) — X]^=o o,jD^ ■ 
As usual [a)j will denote the Pochhammer symbol defined by 

(a)o = 1, (a)j ~ a{a + 1) ■ • • (a + j - 1), for j > 1, a £ C. 

Let /i be a moment functional on the real line, that is, a linear mapping /i : P — M. 
The n-th moment of is defined by /z„ = x"). It is well-known that any moment 
functional on the real line can be represented by integrating with respect to a Borel 
measure (positive or not) on the real line (this representation is not unique [7]). If 
we also denote this measure by /i, we have = J p{x)d^{x) for all polynomial 

p S P. Taking this into account, we will conveniently use along this paper one 
or other terminology (orthogonality with respect to a moment functional or with 
respect to a measure). 

We can associate to the moment functional /x a bilinear form defined in the linear 
space of real polynomials: 

bf] (2.1) B{p,q) = {fi,p{x)q{x)). 

If the moment funcional is represented by a discrete measure supported in a 
countable set X of real numbers, namely fi — X^^g^ /^(^)^^' bilinear form 
above is equal to B{p,q) = J2xex Pi^)1i■'^)^^i^)■ 
Let -ftT be a positive integer. If 



/mo • • • 

defth I (2.2) Qn = det ■ •.. : 

\Mra • • • M2ri 



^0 
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for n = 0, 1, ■ • • , K, we can associate to the moment functional fi a sequence of 
monic orthogonal polynomials (pn) n=o^ Pn of degree n, using the formula 

/ Mo Ml •■• M« \ 



fdpo (2.3) Pn{x) 



1 



■det 



Mn-l 
1 



fJ-n 
X 



M2n 



1 < n < is:, po = 1- 



The polynomials Pm < n < K, satisfy that B{pn,x"^) — 0, ior < m < n and 
B{pn,x'^) ^ 0. Conversely, any polynomial q of degree n, < n < N, satisfying 
B{q,x^) = 0, for < m < n, is equal to p„ up to a multiplicative constant: 
q{x) = apn{x). 

The moment funcional fi can be represented by a positive measure with infinitely 
many points in its support if and only if > 0, for n — 0, 1, • • • , and then the 
associated bilinear form ()2.1|) is an inner product {B{p,p) > for all polynomial 
p^O). 

Favard's Theorem establishes that a sequence {pn)n of polynomials, Pn of degree 
n, is orthogonal with respect to a moment functional if and only if it satisfies a 
three term recurrence relation of the form {p-i = 0) 

XPn{x) = anPn+l{x) + bnPn{x) + C„p„_i(x), n > 0, 

where (a„)„, (&„)„ and (c„)„ are sequences of real numbers with a„c„ ^ 0, n > 1. 
If, in addition, a„c„ > 0, n > 1, then the polynomials {pn)n are orthogonal with 
respect to a moment functional which can be represented by a positive measure, 
and conversely. 



The kind of transformation which consists in multiplying a moment functional fi 
by a polynomial r is called a Christoffel transform. The new moment functional rfi 
is defined by {rfi,p) — {fi, rp). It has a long tradition in the context of orthogonal 
polynomials: it goes back a century and a half ago when E.B. Christoffel (see [5] 
and also studied it for the particular case r{x) = x. 

Given a moment functional and a polynomial r{x) — {x — ai) ■ ■ ■ (x — am), we 
can give, using [36], Th. 2.5, a necessary and sufficient condition for the moment 
functional r/i to have a sequence of orthogonal polynomials. Indeed, write {pn)n=o 
(N finite or infinite) for the orthogonal polynomials with respect to fi and A„ for 
the square matrix 

mataO | (2.4) Ki ^ iPn+iia^))t=i,- ,m.i=o,- ,m-i, n < N -m + 1. 

Then the moment functional r/i has a sequence {qn)n=oi ^ ^ N — m + 1, of 
orthogonal polynomials if and only if det(A„) 7^ for n = 0, • • • ,K. 

For a real number A, the moment functional /i(x + A) is defined in the usual 
way {iJ,{x + X),p) = {n,p{x — A)). Hence, if (p„)„ are orthogonal polynomials with 
respect to /i then {pn{x + A))„ are orthogonal with respect to ^i{x + A). 

We will use the following Lemma to construct orthogonal polynomials with re- 
spect to a measure together with a Dirac delta (for related results see [16], [37]). 

addel Lemma 2.1. From a measure v and a real number A, we define the measure fi by 
H — {x ~ X)^. Assume that we have a sequence {pn)n of orthogonal polynomials 
with respect to ^ and write an = J Pndv. For a given real number M we write 
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hk2 



nip 



def sj 



bet jac3 



p = V + M5\ and define the numbers 



(2.5) 



n> 1 



an-i + Mp„-i{X) ' 

(we implicitly assume that + Afp„_i(A) 7^ 0, n > 1). Then the polynomials 

defined by qg = 1 and qn — Pn + PnPn-i, ore orthogonal with respect to p. 

Proof. We have to prove that J{x — \y qndp = 0, j = 0, • • • , n — 1 and J{x — 
Xrqndp ^ 0. 

From the orthogonahty of {pn)n with respect to p = {x — X)!/ and the definition 
of p we automatically get for j — 1, ■ ■ ■ ,n — 1 

{x - Xyqndp ^ {x ~ Xy~^qndp = 0, 



□ 



1 and 



whatever the numbers n > 1, are. For j = 0, using (|2.5|) . we get 

qndp = J {pn + (3nPn-l)dv + Mpn{X) + /3„M]3„_i(A) 

= an + /3„a„_i + Mpn{X) + /3„Afp„_i(A) = 0. 

Finally, for j = n we have J{x — A)"g„dp — f3n J{x — Xy^^^Pn-idp ^ 0. 

We will need the following technical Lemma. 

Lemma 2.2. Let u be a real number. We consider the polynomials sq.u 
for j >1 

(2.6) s,,u{x)^{-iyll[x + i{u-i)]. 

i=Q 

Then for j > 0, Sj.„(x(a; — u)) = {—x)j{x — u)j. 

Proof. If wc expand {—x)j{x — u)j, we have 

{—x)j{x — u)j = (— l)-'x ■ ■ ■ {x — j + l){x — u) ■ ■ ■ {x — u + j — 1). 

Multiplying the {i + l)-th and the {j + i + l)-th factors in the product above 
{0 < i < j — 1), we get 

(x — i){x — u + i) ^ x{x — u) + i{u — i). 

This gives for {—x)j{x — u)j the formula 

(2.7) {~x)j{x — u)j = {—ly W\x{x — u) + i{u — i)] — Sj;„(a;(x — u)). 



i=0 



3. The method and its main ingredient 



□ 



In this Section, we introduce our method and its main ingredient. Although in 
this paper we will apply it to difference and differential operators and some of the 
classical discrete and classical families of orthogonal polynomials, the method also 
works with g-difference operators and the corresponding g-classical polynomials. 
This is the reason why we introduce the method in an abstract setting. 

The starting point is a sequence of polynomials {pn)n, Pn of degree n, and an 
algebra of operators A acting in the linear space of polynomials. In addition, we 
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assume that the polynomials pn, n > 0, are eigenfunctions of certain operator 
Dp e A. We write (6'„)„ for the corresponding eigenvalues, so that Dp{pn) = 6'„p„, 
n > 0. 

Given a sequence of numbers (/3n)n, we define a new sequence of polynomials 
iqn)n by go = 1 and 

qnCC I (3.1) qn = Pn + PnPn-1, n>l. 

Assuming certain structure for the sequence (/3n)n, our procedure provides a new 
operator Dq G A for which the polynomials ((?„)„ are eigenfunctions. The structure 
for the sequence (/3n)„ will be determined by the main concept involved in our 
method: 2?-operators. 

We have two different types of P-operators, depending on whether the sequence 
of eigenvalues (0„)„ is linear in n or not. 

Given a sequence of numbers (e„)„, a I?-operator of type 1 associated to the 
algebra A and the sequence of polynomials is defined as follows. We first 

consider the operator C : P — P defined by linearity from C,{pn) — £nPn-i- Using 
that dgr(C(p)) < dgr(p) for all polynomial p, we next define the operator V by 

P : P ^ P 



defTo 



(3.2) V{p)^Y.^-iy^^C'{p). 

We then say that D is a P-operator if P e ^. 

In the next Sections, we will display a P-operator of type 1 for Charlier and 
Laguerre polynomials and two 2?-opcrators of type 1 for Meixner and Krawtchouk 
polynomials (see Lemmas 14.11 15.11 16.11 and 18. ip . 

Once we have introduced the main concept of P-operator of type 1 defined by the 
sequence of numbers (e„)„, we are ready to go on with our method. The estructure 
mentioned above for the sequence (/3„)„ (which defines the polynomials ((?„)„ p.ip ) 

is then /3„ — e„ — ^ " \ i where P is any arbitrary polynomial. The details of 

P[On-l) 

the method are included in the following Lemma. The Lemma is an elementary 
but powerful result: its proof is just an easy computation, but it provides higher 
order differential, difference or g-difference operators for a large class of polynomials 
(some of which turn out to be orthogonal with respect to a measure) . For the case 
of differential operators, these polynomials include the Krall-Laguerre-Koornwinder 
families orthogonal with respect to the measures (I1.17p . 

fllv Lemma 3.1. Let A and {pn)n be, respectively, an algebra of operators acting in 
the linear space of polynomials, and a sequence of polynomials, pn of degree n. We 
assume that {pn)n are eigenfunctions of an operator Dp € A, that is, there exist 
numbers On, n > 0, such that Dp{pn) = dnPn, n > 0. We also have a sequence of 
numbers (e„)„ which defines a V-operator of type 1 for {pn)n and A (see 13.S\} ). 
For an arbitrary polynomial P2 such that P2{9n) ^ 0, n > 0, we define the sequences 
of numbers (7„)„ and (A„)„ by 

(3.3) 7„+i = P2{9n), n>0, 
def jg I (3.4) A„ - A„_i = 7„, n > 1, 
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and assume that there exists a polynomial Pi such that A„ = Pi(0„). We finally 
define the sequence of polynomials {qn)n by = 1 and 



defqng (3.5) g„ = p„ + /3„p„_i, n > 1, 

where the numbers n> 0, are given by 



In+l 



defbetng (3.6) /3„ — e„ , n > 1. 

1 7« 



defD 



Then Dq{qn) = Xnqn where the operator Dq is defined by 
(3.7) Dq=Pi{Dp)+VP2iDp). 
Moreover Dq £ A. 

Proof. Since Dp{pn) = 0„Pn, we have that 

Pl{Dp){pn) = Pl{On)Pn = A„p„, 
P2{Dp)(pn) = P2{0n)Pn = In+lPn- 

The definition of Dq (|3.7p then gives 

DqiPn) = A„p„ + -1n+l'D{pn)- 

Using the definition of the 2?-operator V (|3.2p . we have 



Dqipn) - A„p„ + 7„+i ^(-l)^'+iC'(Pn) 

n 

= A„p„ + ^7„+i(-l)-' + ^£„ • ■ • En-j + lPn-j- 
J = l 

Taking into account the definition of the polynomials ((?„)„ p.Sp . we get 

Dq{qn) = Dq{pn) + PnDq{pn-l) 

n 

= A„_p„ +7n+ienPn-l + ^ 7„+i (- 1 )-'+^ £„ • • • e„-j+lP„-j 

n-1 

+ /3„A„_ip„_i + /3n ^ 7„(-l)-'+^e„_i • • • e„_jp„_i_j- 



eql (3.8) = XnPn + {PnXi-1 + 7ri+ien)Pn-l + ^ an,jPn-j, 

where 

a„j = (-l)^+^(7„+ie„ • --en-j+i - /3„7„e„_i • --en-j+i), 2 < j < n. 
Using the definition of A„ and /3„ f p.4p and p.6p ). we easily have 

A„/3 

a„.j =0, 2 < j < n. 

From where (|3^ gives that Dq{q„) = XnPn + XnPnPn-i = Xnqn- 

Since 2? is a P-operator for the polynomials {pn)n and the algebra A, we have 
that V € A, and then D„ G ^ as well. □ 
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Actually, a slight modification of the Lemma provides a huge family of operators 
in A for which the polynomials {qn)n are eigenfunctions. Indeed, consider any poly- 
nomial G for which there exists a polynomial Pi satisfying that Pi(0„) — Pi(0„_i) = 
G{0n-i)P2{dn-i) (when the eigenvalues On are linear in n, that polynomial Pi al- 
ways exists for each polynomial G). Proceeding as in the previous Lemma, it is 
easy to see that the operator Dq c defined by Dq,c — Pi{Dp) + G{Dp)'DP2{Dp) 
belongs to A and satisfies Dg^ciQn) = Pi{On)qn- 



defTo2 



fl2v 



def qng2 



defbetng2 



We will apply the version of our method established in Lemma 13.11 when the 
eigenvalues {9n)n of the polynomials {pn)n with respect to the operator Dp are 
linear in n. It is easy to see that if {On)n are, say, quadratic in n, the hypothesis of 
Lemma l3. II can not fulfil. Indeed, if jn+i is a polynomial in 0„ and 6n are quadratic 
in n, then 7„ is a polynomial in n of even degree, and the sequence A„, defined 
by A„ — A„_i = 7„, is a polynomial in n of odd degree. Hence, it can not be a 
polynomial in the quadratic sequence 9n- 

To fix this situation we need to introduce a modified version of the 2?-operator 
defined in (|3.2|) . Two sequences of numbers (e„)„ and (cr„)„ are now the starting 
point. A P-operator of type 2 associated to the algebra A and the sequence of 
polynomials {pn)n is defined as follows. As before we first consider the operator 
C : P — ^ P defined by linearity from ({pn) = £nPn-i- We next define the operator 
2? again by linearity from 



(3.9) 



V{pn 



1 



lP„ + > (-l)^ + V„+i_,C^(Pn) 



We then say that P is a P-operator liV £ A. 

When the sequence ((T„)„ is constant, a D-operator of type 2 reduces to a 2?- 
operator of type 1 (up to constants). 

In Lemmas 17.21 and 18. 7[ we display four D-operators of type 2 for the Hahn 
polynomials and two 2?-operators of type 2 for the Jacobi polynomials, respectively. 

We can now establish the version of our method for 2?-operators of type 2. 

Lemma 3.2. Let A and {pn)n be, respectively, an algebra of operators acting in 
the linear space of polynomials, and a sequence of polynomials, pn of degree n. We 
assume that {pn)n are eigenfunctions of an operator Dp g A, that is, there exist 
numbers On, n>0, such that Dp{pn) = 9nPn, n>Q. We also have two sequences 
of numbers (£„)„ and {crn)n which define a D-operator of type 2 for {pn)n and A 
(see \3.9\) ). For an arbitrary polynomial P2 such that P2{9n) ^ Q, n>Q, we define 
the sequences of numbers {'jn)n and (A„)„ by 



(3.10) 
(3.11) 



A. 



7«+i 

Ari-l 



P2{9n), 



n > 0, 
n > 1, 



and assume that there exists a polynomial Pi such that Xn+i + A,, 
finally define the sequence of polynomials {qn)n by Qo = ^ and 



Pi{9n)- We 



(3.12) 

where the numbers /3„ 
(3.13) 



qn = Pn + PnPn- 

n>Q, are given by 

In+l 



1, n>l, 



n>l. 



In 
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defD2 



chw 



repnch 



sodech 



sdch 



ITch 



Then Dq{qn) — Xnqn where the operator Dq is defined by 



(3.14) 



Moreover Dq G A. 



Proof. The Lemma can be proved proceeding in the same way as in the proof of 
Lemma 13.11 □ 



We remark that if an operator D is a P-operator for the polynomials (p„)„ and 
the algebra A then it is also a 27-operator for the sequence of polynomials Pn — OnPn, 
whatever the real numbers a„ 7^ 0, n > 0, are. Indeed, it is straightforward to see 
that if T) is defined from {pn)n by the sequence {Sn)n then V is the defined from 
{P7i)n by the sequence 



defchp (4.1 



Using this fact, we can associate P-operators to a measure having orthogonal poly- 
nomials. Indeed, we say that an operator I? is a I?-operator for a measure fj, if it is 
a P-operator for a sequence {pn)n of orthogonal polynomials with respect to fj,. We 
point out that the sequence of numbers (£n)n, which defines the ©-operator from 
the sequence of polynomials (pn)n, depends on this sequence of polynomials. In 
other words, given a measure /it the corresponding set of P-operators only depends 
on the measure and it is independent of the orthogonal polynomials normaliza- 
tion, but the associated sequences of numbers (e„)„ depend on this normalization. 

4. ChARLIER CASE 

For a 7^ 0, we write (c"j)„ for the sequence of Charlier polynomials defined by 

1 " / \ / \ 



3 J \J 



J'- 



(that and the next formulas can be found in [4j, pp. 170-1; see also [13]j PP) 247-9 or 
[35] , ch. 2). The Charlier polynomials are orthogonal with respect to the measure 



(4.2) 



Pa 



a ^ 0. 



The measure pa is positive only when a > 0. 

The three term recurrence formula for (cj'jn is {cti = 0) 

(4.3) xcn — {n + l)cn+i + {n + a)cn + acn-i, n>0 

(to simplify the notation we remove some parameters in some formulas) . They are 
eigenfunctions of the following second order difference operator 

(4.4) Da^xS-i-{x + a)So + aSi, !)„«) = -«<, n > 0. 
They also satisfy the simple difference relation 

(4.5) A(c„)=c„_i. 

We first identify a I?-operator for the Charlier polynomials. 

Lemma 4.1. For a ^ 0, the operator V defined by 13. S\) from the sequence = 1, 
rt > 0, is a "D-operator for the Charlier polynomials (cj^)n {4-1^ and the algebra A 
il.lO\) of difference operators with polynomial coefficients. More precisely I? = V. 
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Proof. Indeed, since A(c^) = C^-i (|4.5p . we have that C = A. Taking into account 
that AQ) = (J_J,m>l, we find 

\mj ^-^ \m I ^-^ \TO— 7 

\ rix , \ / \ / 1 



This gives P = V G yl. 



□ 



Lenima l3.1l and the 2?-operator for the Charlier polynomials found in the previous 
Lemma, automatically produce the large clase of polynomials satisfying higher order 
difference equations displayed in Theorem 11.11 in the Introduction. We are now 
ready to prove it. 

Proof. Indeed, the theorem is a straightforward consequence of the Lemma [01 We 
have just to identify who the main characters are in this example. 

If we write e„ = 1, Lemma |4 . 1 1 gives that the sequence defines a P-operator 

for the Charlier polynomials (c°)„ and that P = V. 

The eigenvalues of (c"j)„ with respect to Da arc 6'„ = — n, hence the definition of 
(7n)n, (An)n and the relationships between the polynomials P\ and P2 automati- 
cally give 

An — A„_i = 7„. 

The Lemma [3. II now shows that the polynomials g„, n > 0, defined by (|1.14p are 
eigenfunctions of the difference operator D = Pi{Da) + '^P-ziDa) with eigenvalues 

An- 

Since Pi and P2 are polynomials of degrees fc + 1 and fc, respectively, and Da 
has genre (—1, 1), the operators D{Pi) and D{P2) have order 2k + 2 and 2k and 
genre (— fc — l,fc + 1) and {—k,k), respectively. Write now /_i = x and fi = a 
for the coefficients of 5_i and Si in Da, and ui, U2 for the leading coefficients of 
the polynomials Pi and P2 , respectively. A straightforward computation gives that 
the coefficients of 5-k-i and Sk+i in the operator D are f^i{x — 1) • • • f-i{x — 
k){uif^i{x) — U2) and uifi{x) ■ ■ ■ fi{x + fc), respectively. We can see that both 
are different to zero. This shows that the operator D has order 2fc + 2 and genre 
(-fc- l,fc+l). 

□ 

We next find some particular choices of P2 for which the sequence of polynomials 
{Qn)n (|l-14p is orthogonal with respect to a measure. 

For k a positive integer, let pk^a be the measure defined by (a 0) 

(4.6) pk,a ^{X+1)---(X + k)pa{x + fc + 1) 

°° x+k+1 

^(-im-.-.+i: ^^^^ ..^ 

a:;— 

To simplify the notation we will sometimes remove the dependence of a and fc and 
write p = pk,a- 

In order to find the orthogonal polynomials with respect to pk,a we need the 
following Lemma. 
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chx+k+1 



chxx Lemma 4.2. For a positive integer k and a non-null real number a, we have 
(4.7) (p,,„<) = (-ire"fc!c,-"(-n-l). 

Proof. Write 

^« = (pfe,a,0, n > 0, 
C„ = (-irfc!e"c-'^(-n-l), n>0. 
The three term recurrence relation (I4.3P for (c^)„ gives that 



meccla 



meccl 



mecc2a 



mecc2 



ei + (a + fc + l)eo-a 



fc+1. 



0, 



(4.8) 

(4.9) (n + l)e„+i + (n + a + fc +!)$„ + aC„-i 0, n > 1. 
We now prove that also 

(4.10) Ci + (a + ^- + l)Co-a''+'e'^ = 0, 

(4.11) {n + l)Cn+i + {n + a + k+l)Cn + aCn-i = 0, n>l. 
Indeed, for n > 1 we have 

{n + l)C„+i + ia + n + k + 1)C„ + - (-l)"+ifc!e" 

X ((n + l)cj:"(-n - 2) - (71 + a + fc + l)cj:''(-n - 1) + ac-"(-n)) . 

Then, by writing x = — n — 1 in the second order difference equation (|4.4p for the 
Charlier polynomials (c^")fc, we find 

(n + l)C„+i + {a + n + k + 1)C„ + aCn-i = 0. 

Proceeding in a similar way for n — 0, and using that c^^iO) = a'', we get (|4T0l) . 

Since the sequences (^n)ri and (Cn)n satisfy the same recurrence relation, it is 
enough to prove that = Co- Write then rjo = {pa{x + 1) , I) = e°- , and 

Vk - {Pa{x + k + l),{x + l)---{x + k)), k>l. 

From the definition of pk.a we have rjk ~ Co- Using that xpa{x) — apa{x — 1), we 
get 

Vk = {pa{x + k + 1),{X + 1) ■ ■ ■ {X + k)) ^ {pa{x + k),x{x + I) ■ ■ ■ {X + k - I)) 

= {apa{x + k - l),{x + I) ■ ■ ■ [x + k - I)) - k{pa{x + k),{x + l)---{x + k- 1)) 

= a^{pa, 1) - kT]k-i = a'^e" - krjk-i. 

This means that the sequence rjk, k > 1, is characterized by the recurrence relation 

rjk + krjk-i = a^e", k>l, 

with initial condition rjo — . The proof will finish if we prove that also the 
sequence [e°'k\c~^°'{—V))h satisfies the same recursion. But, by applying k — \ times 
the three term recurrence relation for the Charlier polynomials (c~^°'^k (see ()4.3p ). 
we get 



Cr(-l)+C^-l(-l) = ^(Cfe-l(-l)+Cfe^2(-l)) 



„fc 1 „fc 

^(cr(-l) + c-(-l))^^. 

□ 



We are now ready to prove that the orthogonal polynomials with respect to pk^a 
4.6p are a particular case of Theorem 11.11 
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lm51 



eigchk 



def delch 



def betch 



def qch 



defplch 
defp2cli 



jodch 



Theorem 4.3. For k > I, let a be a non-null real number satisfying that ""{—n) ^ 
0, n > 1, where c^,"" , k > 1, are Charlier polynomials (see ^.1^ ). We define the 
sequences of numbers (A„)„, (7„)„>i and (/3„)„>i by 



(4.12) 
(4.13) 
(4.14) 



ln+1 



n > 1, 



In 



and the sequence of polynomials {qn)n by qo = 1, and 

(4.15) qn{x) = <(a;) + /3„<_i(a:), n > 1. 

Then the polynomials (q„)„ are orthogonal with respect to the measure p {4^.6^ . 
Moreover, consider the difference operator of order 2k + 2 and genre {—k — 1, fc + 1) 
defined by 

D = Pi{Da)+VP2{Da), 

where Da is the second order difference operator for the Charlier polynomials {4-4^ 
and Pi, P2 are the polynomials of degrees k + 1 and k, respectively, defined by 

(4.16) Pi{x) ^ -c^^,{x), 

(4.17) P2{x)=c^^ix-1). 
Then D{q„) = A„g„. 

Proof. In the notation of Lemma |2.H we have A = — fc — 1, ly = Pk,a, M = of^^^ Pai 
Pn = c° and M = 0. Lemma l42l gives that a„ = {'^,Pn) = (-l)"e''fc!c^°(-n - 1). 
The orthogonahty of the polynomials ((/„)„ with respect to pk.a is now a conse- 
quence of Lemma [O and (|4.13p . (|4.14p . 

The second part of the Theorem is a straightforward consequence of the Theorem 
O We have just to check that P2{x) ^ Pi{x - I) - Pi{x): 



Pi{x~l)-Pi{x) 



-cll,{x-l)^A{cll,{x-l)) 



cr(^-l) 



P2{x) 



□ 



Except for an appropriate choice of the polynomial P2 (as before), the polyno- 
mials {qn)n (|1.14p in Theorem 11.11 seem not to be orthogonal with respect to a 
measure. Anyway, in a subsequent paper |10| . we prove that they enjoy certain 
orthogonality property which implies that they satisfy a higher order recurrence 
relation of the form 

ko + l 

{x + l){x + 2) ■ ■ ■ [x + ko + l)qn{x) = ^ anjqn+j{x), 

j = -ka-l 

where fcp < deg(P2)- 

Conjecture 2 in the Introduction for the Charlier weight pa and the finite set 
F = {1, 2, • • ■ , fc} can be deduced as a Corollary of the previous Theorem. 

Corollary 4.4. For a positive integer k, write F = {1,2,- •• , fc} and Pf{x) — 
rijG-F^^"-^)" ^•'SMme that a ^ and that c'j^°'{—n) 7^ 0, n > 1. Then, the orthogonal 
polynomials with respect to the measure pppa are common eigenf unctions of a {2k -\- 
2)-order difference operator of genre (— fc — l,fc -|- 1), where the measure pa is the 
Charlier measure li4-0)- 
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remch 



condepch (4.18) 



Proof. The proof is an easy consequence of the previous Theorem, taking into 
account that the measure pppa is equal to the measure p shifted by A; + 1: pppa = 
p{x-k-l). □ 

Remark 4.5. The most interesting case in the previous Theorems is when the 
new measures pk,a (|4.6p in Theorem 14.31 and ppPa in CoroUarv 14.41 are positive. 
This obviously happens when a > and k is even. In that case, the hypothesis 
c^°(— n) 7^ 0, n > 1, in Theorem 14.31 and CoroUarv 14.41 alwavs fulfil. Indeed, by 
definition of the measure pk.a (|4.6p we have that 

PkM = (a; + 1) • • • (a; + k)pa{x + k + l) 

where pa{x + k + 1) is the Charlier weight (|4.2p shifted by —A; — 1. 

That means that the measure pk.a is a Christoffel transform of the measure 
Pa{x + k + 1). Write A„, n > 0, for the matrices 



— (c„+j_i(*))ij-=l 



Hence, according to (j2.4p in the Preliminaries, the measure pk.a has a sequence 
('7n)n of orthogonal polynomials if and only if det A„ ^ 0, n > 0. Now for a > 
and k even, the measure pk.a is positive with infinitely many points in its support. 
So it always has a sequence of orthogonal polynomials, and then det A„ ^ 0, n > 0. 
In [9] we have found a very nice closed-form expression for the determinant of the 
matrices A„ (|4.18p : 



chconjdet (4.19) 



det(A„) = 



(n + j-1)! 



cr(-")- 



We can then conclude that the hypothesis c^°'{—n) 7^ 0, n > 1, in Theorem 
and CoroUarv 14.41 alwavs fulfil for a > and k even. 

There is a similar formula to (|4.19l) for each of the examples studied in Sections 
5, 6 and 7 of this paper. Actually, they seem to be particular cases of a surpris- 
ing symmetry for certain Casorati determinants associated to the classical discrete 
orthogonal polynomials (see [9]). 



defmep 



sodeme 



trme 



5. MeIXNER CASE 

For a 7^ 0, 1 we write {ra'^'^)n for the sequence of Meixner polynomials defined 



by 

(5.1) 



mT{x) = (-1) 



n \ ^ —7 



(we have taken a slightly different normalization from the one used in j4], pp. 175- 
7, from where the next formulas can be easily derived; see also [23], pp, 234-7 or 
[35) . ch. 2). Meixner polynomials are eigenfunctions of the following second order 
difference operator 
(5.2) 

Da.c = x5-i-[{l + a)x + ac\5a + a{x + c)5i, Da^drn'^f) n(a-l)m°''', n > 0. 
When a 7^ 0, 1, they satisfy the following three term recurrence formula (m_i = 0) 



(5.3) 



xrrLr, 



a{n + 1) 



(a + l)n + ac 



1 



-m„_i, n > 
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(to simplify the notation we remove the parameters in some formulas). Hence, for 
a 7^ 0, 1 and c^O,— 1,— 2,---, they are always orthogonal with respect to a moment 
functional Pa.c, which we normalize by taking {pa.c, 1) = r(c). For < |a| < 1 and 
c / 0, —1, —2, ■ • • , we have 

[miF| (5.4) pa^,, = (1 - af 2^ ^— '-5^. 

The moment functional pa,c can be represented by a positive measure only when 
< a < 1 and c> 0. 

Meixner polynomials satisfy the following identities 

S (5-5) A(mr) = 



sdm2 



sdm2b 



(5.6) mr+H-^-i) = E«'-''<A(^)' 



J=0 



(5.7) mT+'{x)=Y,m:'l^{a 



(5.8) ^m:;l,ix)^mr{x)-m^^'^-'ix + l). 

For Meixner polynomials we have found two different 2?-operators. 

ITme Lemma 5.1. For a 0,1, consider the Meixner polynomials (mj^''^)n i5.1\) . Then, 
the operators Vi, i — 1,2, defined by iS.S^) from the sequences (n>0) 

vemel (5.9) £„,1 = —1, 

veme2 (5.10) £„,2 = — 

are V-operators for {m'^''^)n and the algebra A 11.10\) of difference operators with 
polynomial coefficients. More precisely 

domel I (5.11) 2?i = — ^A, 

' ' 1 — a 

dome2 I (5.12) X'2 = ^ ~ 



1-a 

Proof. Indeed, the sequence (|5.9[) defines the 2?-operator p. 21) 



i=i j=i i=o 

Using (1221), (lOI) and (1531) . we get 

Pi(m:^^=) = ™:5'^(x) - m«'^+i(x) = m'^fix) - im:f+^(x) - m'},^'^{x + 1) 

= T^A(™:^'-). 
1 — a 

The proof for 2?2 is similar and it is omitted. 

□ 

Each one of the P-operators displayed in the previous Lemma together with 
Lemma 13.11 gives rise to the corresponding class of polynomials satisfying higher 
order difference equations. 
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th5 . Imel 



beme 



qnme 



Theorem 5.2. Let Pi he a polynomial of degree k + I, k > 1, and write P2{x) = 
Pi{x + a — 1) — Pi{x) (so that P2 is a polynomial of degree k). We assume that 
P2{n{a — 1)) 0, n > 0, and define the sequences of numbers 



(5.13) 
(5.14) 



'fn+i = P2{n{a - 1)), n>0, 
A„ = Pi(n(a-l)), n>0, 

Pn.i = n>l, 

In 

o 1 7«+l ^ -, 

Pn.2 = n>l, 

a 7« 



and the sequence of polynomials go.i ~ Ij o'^^ for n > 1 
(5.15) qn,iix) = m'^^''{x) + l3n,im'^^liix), i 



1,2, 



where m"^''^ is the Meixner polynomial defined by 15.1]} (a ^ 0, 1 Consider the 
second order difference operator Da,c i5.2\) with respect to which the Meixner poly- 
nomials are eigenf unctions. Write finally Di, i ~ 1,2, for the difference operators 
of order 2k + 2 and genre {—k — 1, /c + 1) 



Di=Pi{Da,c) 



-AP2iDa,. 



1 - a 

D2 = Pl{Da^c) + -^VP2(-Da,, 



1 



Then Di{qn,i) = A„g„,i, n > 0, i ^ 1,2. 



Proof. The theorem can be proved as Theorem 11.11 using Lemma 13.11 and the 2?- 
operators in Lemma l5. II 

□ 



tmew2 



tmewb2 



Imelx 



melx+k+1 



5.1. Meixner I. For a convenient choice of the polynomial Pi, the polynomials 
(9n,i)n in the previous Theorem turn out to be orthogonal with respect to a moment 
functional. Indeed, for k a positive integer, a ^ 0,1 and c k + 1, k, , k — 1, ■ ■ ■ , we 
consider the moment functional 



(5.16) 



PkM.c = {x + C - 1) ■ ■ ■ {x + C - k)pa,c-k-l. 



For < |a| < 1 and c ^ k + l,k, - ■ ■ the moment funcional pk,a,c can be represented 
by the measure 



(5.17) 



Pk.a.c = (1 - a) 



a^r(a; + c) 



^ (x + c — k — l)x\ 

x—Q ^ ' 



5x- 



To simplify the notation we will sometimes remove the dependence of a, c and k 
and write p = pk,a,c- 

The following Lemma will be needed to find the orthogonal polynomials with 
respect to pk,a,c. 



Lemma 5.3. For a 0,1, c 7^ fc + 1, fc, • 

{-l)''kir{c^ 



(5.18) 



• and n > 0, we have 

k — l)m]/°'' '^^^(— n — 1) 



{1-ar 
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Proof. We proceed in an analogous way to Lemma l4?2] Write 

^n = {Pk,a,c,ml'''), n>0, ^-1 = "'"*''^ 



Cn 



(-l)'=fc!r(c - fc - l)m^/"'""+^(-n - 1) 



n > -1. 



(1-a)^ 

The three term recurrence relation ()5.3p for (?7i^j''^)„ gives that 

me erne 1 | (5.19) a„,^„+i + (6„ + c - /c - 1)^„ + c„^„_i = 0, n > 1, 

where a„, and c„ are the recurrence coefficients of the Meixner polynomials 
(m°'°)„. Using that {x + c)pa,c — Pa,c+i/(l — a), the normalization for the Meixner 
weight and the value of we also have (j5.19p for n = 0. 
We now prove that also 

niecme2 | (5.20) a„C„+i + (6„ + c - /c - l)Cn + c„C„_i = 0, n > 0. 

Taking into account the definition of C„, this is equivalent to 

a{n + l)m^/"'~"+'(-n - 2) - ((a + l)(ji + 1) - k{l - a) + c - 2))m^/"'""+'(-7i - 1) 

+ (n + c- l)TO^/°'"'=+^(-?i) = 0, n>0. 

But this follows straightforwardly by writing x — ~n — \ in the second order 
difference equation for the Meixner polynomials ("^J,^"' '^^^)fe (see (|5.2p ). 

Since the sequences {£,n)n and (Cn)n satisfy the same recurrence relation ( (|5.19p . 
(|5.20p ). it is enough to prove that = C-i and = Co- 
A straightforward computation from (|5.ip shows that = 
In order to prove that also i^o — C,o^ we proceed as follows. Write 

rik,c = {pa,c-k-i, (x + c - 1) • • • (a; + c - k)). 

From the definition of pk,a,c we have rik,c — Co- Using that {x + c)pa.c — Pa,c+i/(l — 
a), we get 

Vk.c = (Pax^fc-i, (x + c-l)- --{x + c- k)) 
{Pa,c-k: ix + c-2)-- - (x + c-k)) 



1 - a 

-J- + krik-i,c-i = fc»7fc-i,c-i + 



+ k{pa,c-k-l, (X + C - 2) ■ 

r(c-i) 



■(x + c-k)) 



(1-a)* 



(l-a) 

On the other hand, if we write 

(-l)'=fc!r(c - fc - l)m^/"'"'+^(-l) 



Tk,. 



Co 



(1 -fl)*^ 



we have using (15.51) and (15.11) 

(-l)'=fc!r(c- fc 



Tk,, 



kTk-l,c-l = 



(l-a)'= 
(~l)'=fc!r(c- fc 



-(m 



l/a,-c+2 



(-l) + (l-a)mi./j'r+^(-l)) 



i),„;/..-,„,.Qi_;), 



Proceeding by induction on k it is now easy to prove that rjk c — Tk c, and then 
Co = Co. 

□ 
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lm62 



def delme2 



defbetme2 



We are now ready to prove that the orthogonal polynomials with respect to pk,a,i 
()5.16|) are a particular case of Theorem 15.21 



Theorem 5.4. For k > 1, let a and c be real numbers satisfying that a ^ 0,1, 
c ^ A; + 1, fc, fc — 1, • • • and rn^/'^' '^^'^{—n) ^ 0, n > 1, where rn]/'^' '^^'^ , k > \, 
are Meixner polynomials (see \5.1\) ). We define the sequences of numbers (A„)„, 
(7n)n>i and (/3„)„>i by 



a- 1 



(5.21) 
(5.22) 



l/a,-c+2 



l/a,-c+l, N 



7«+i 



and the sequence of polynomials (<?„)„ by go = 1, o,nd 
def qme2 I (5.23) (7„ (a;) = m°^'=(x) + /3„m";:!;^ (x), n > 1. 



Then the polynomials {qn)n are orthogonal with respect to the moment functional 
Pk,a,c \5.16]) . Moreover, consider the difference operator of order 2k + 2 and genre 
{—k — 1, fc + 1) defined by 



D = Pi{Da,c) 



l-a 



where Da^c is the second order difference operator for the Meixner polynomials 
and Pi and P2 are the polynomials of degrees k + I and k, respectively, defined by 



Pi{x) = 



1 



a- 1 



l/a,-c+l 



X 



P2{x) = rUf. 



l/a,-c+2 



a-1 



a- 1 
- 1 



Then D(g„) = A„g„. 

Proof. The Theorem can be proved in an analogous way to Theorem 



□ 



lm623 



Conjecture 1 in the Introduction for the Meixner weight pa,c and the finite set 
F = {1, 2, • • • ,k} can be deduced as a Corollary of the previous Theorem: 

Corollary 5.5. For a positive integer k, write F — {1, 2, • • • , fc} and P^c-f{x) = 
Y\j^F^^ + c + j). Assume that a ^ 0, 1, c ^ 0, —1, —2, • • • and that 



l/a, — c— fc+l/ \ / n \ 1 



m 



Then the orthogonal polynomials with respect to the measure p -c-F P a, c are common 
eigenf unctions of a {2k -\- 2.) -order difference operator of genre {—k— 1, fc+ 1), where 
Pa^c is the moment functional for the Meixner polynomials. 

Proof. The proof is an easy consequence of the previous Theorem, taking into 
account that the measure p-c-FPa,c is equal to the measure Pk.a.c+k+i- D 
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tmew 



tmewb 



lm61 



eigmek 



def delme 



def betme 



def qme 



defplme 



defp2me 



5.2. Meixner II. There is other choice of the polynomial Pi for which the poly- 
nomials {qn,2)n in Theorem [52] are also orthogonal with respect to a moment func- 
tional. Indeed, for k a positive integer, a ^ 0, 1 and c ^ fc + 1, fc, • • ■ let pk.a.c be 
the moment functional defined by 



(5.24) 



Pk,a,c = (a; + 1) • • ■ (x + k)pa^c-k-i{x + fc + 1). 



In particular, for < |a| < 1 and c ^ k + l,k, ■ • ■ the moment funcional pk,a,c can 

be represented by the measure 

(5.25) 



Pk,a,c = (1 - a) 



;-fc-l 



(-i)^fc!r(c - k - i)5^k-i + f: 



x=0 



l)x\ 



To simplify the notation we remove the dependence of q, c and k and write p — Pk^a^c 
(the proofs of the following theorems are similar to the ones in the previous Sections 
and are omitted). 



Lemma 5.6. For a ^ 0,1, c ^ k + l,k, ■ ■ 

(-l)'=A:!r(c 



{pk,a,c,m'^i' 



and n > 0, we have 
k-l)ml'~'+^{-n-l) 



iTi — k 



Theorem 5.7. For k > 1, let a and c be real numbers satisfying that a ^ 0,1, 
c ^ fc + 1, fc, , fc — 1, • • • and m^'^'^^^(— n) 7^ 0, n > 1, where to^'~'^^^, k > 1, 
are Meixner polynomials (see i5.1\) ). We define the sequences of numbers (A„)„, 
(7n)n>i and (/?„)„> 1 by 



(5.26) 
(5.27) 
(5.28) 



An — — 



a a,-c+l 
-771 



a- 1 
a.-c+2 

1 7n+l 
a 7n 



fc+1 



arid i/ie sequence of polynomials {qn)n by qa = 1, and 

(5.29) qn{x) = m.'if{x) + /3„7n^;f ^(x), ri > 1. 

T/ieTi the polynomials {qn)n are orthogonal with respect to the moment functional 
p \5.24^ . Moreover, consider the difference operator of order 2k + 2 and genre 
(— fc — 1, fc -I- 1) defined by 

1 



D = Pi{Da,c) 



-VP2(i?a,c), 



1 -a 

where Da.c is the second order difference operator for the Meixner polynomials 115. 2(1 
and Pi and P2 are the polynomials of degrees k + I and k, respectively, defined by 

X 



(5.30) 
(5.31) 

Then D{q„) = A„g„ 



Pi{x) = - 



a- 1 



a.-c+l 
^k+1 



P2{X) 



a, — c4-2 



a - 1 
- 1 



Conjecture 2 in the Introduction for the Meixner weight po,c and the finite set 
F = {1, 2, • • • , fc} can be deduced as a Corollary of the previous Theorem: 



P-OPERATORS, ORTHOGONAL POLYNOMIALS AND DIFFERENCE EQUATIONS 23 



Corollary 5.8. For a positive integer k, write F — {1,2,- •• , fc} and Pf{x) = 
Yijepi^^j)- Assume that a ^ 0,1, c ^ 0, —1, ~2, ■■ ■ and that 

rri^ '^^^^^ {—n) ^0, n>l. 

Then, the orthogonal polynomials with respect to the measure pFPa,c are common 
eigenf unctions of a (2k -\- 2) -order difference operator of genre (— fc— l,fc+l), where 
Pa,c is the moment functional for the Meixner polynomials (rn°''^)„. 

Remark 5.9. Let us note that in some cases the measure pFPa,c can be positive 
even though the measure pa.c is not. Indeed, that it the case, for instance, when k 
is odd, < a < 1 and —2j — 1 < c < — 2j, j ~ 0, - ■ ■ ,{k — l)/2, or when k is even, 
< a < 1 and -2j < c < -2j - 1, j = 1, • • • , k/2. 

6. KrAWTCHOUK CASE 

For a ^ 0, —1, we write (fcJJ'^)ri for the sequence of Krawtchouk polynomials 
defined by 

(6.1) fc-^(.) ^ i-y(-l)"+.-(^ + °)'''" i-r^U--UN-n)^.,_ 

j=o 

Krawtchouk polynomials are eigenfunctions of the following second order difference 

operator (n > 0) 

(6.2) 

Da.N ^ xS-i-{x-a{x~N+l))Sfi'a{x-N+l)5i, L>a,Ar(C^) = ~n(l+a)k^'^ . 

For a 7^ 0, — 1 and N ^ 1,2,---, they are always orthogonal with respect to a 
moment functional Pa,N, which we normalize by taking {pa.N, 1) = 1- When N is 
a positive integer and a > 0, the first N polynomials are orthogonal with respect 
to the positive Krawtchouk measure 

, , r(iv) 

^'^■"^ = (l + a)^-i TiN-x)xl^^- 

The structural formulas for (fcj^'^)n can be found in [3S], pp. 30-53 (see also [23] . 
pp, 204-211). 

Proceeding as in the previous Sections, we have identified two I?-operators for 
Krawtchouk polynomials (the proof is omitted). 

Lemma 6.1. For a ^ 0, —1 and G M, consider the Krawtchouk polynomials 
(fc°'^)„ id. Ip . Then the operators T>i, i — 1,2, defined by I13.2\) from the sequences 
(n > ) 

£ns = 1/(1 + a), 
£n,2 = -a/(l + a), 

are T>-operators for {k!^'^)n and the algebra A hi. 10]) of difference operators with 
polynomial coefficients. More precisely 

1 + a 

1 + a 
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th5 . Ikr 



bekr 



qnkr 



tkrw 



tkr2 



Each one of the I?-operators displayed in the previous Lemma together with 
Lemma 13.11 gives rise to the corresponding class of polynomials satisfying higher 
order difference equations. However, we only proceed with the first 2?-operator 
because due to the symmetries of Krawtchouk polynomials, the examples generated 
from the second D-operator are essentially the same: if we write Qn^j^, InliN^ 
respectively, for the polynomials defined by (j6.5|) below from the first and the second 
©-operators above, it is not difficult to see (using that fc°'^(x) = { — l)"kn^°''^{N — 
1 — a;)) that 

'?i!U(^) = (-i)>i'UA'(-^+^-i) 

(using the first 2?-operator we will prove Conjecture 2 in the Introduction for the 
Krawtchouk polynomials, while using the second P-operator one can similarly prove 
the equivalent Conjecture 2' in 8 for the Krawtchouk polynomials). 

Theorem 6.2. Let Pi be a polynomial of degree k + I, k > 1, and write P2{x) = 
Pi (x — 1 — a) — Pi{x) (so that P2 is a polynomial of degree k). We assume that 
^2(^'T-(1 + a)) 7^ 0, n > 0, and define the sequences of numbers 

7„+i = P2(-^^(l + a)), > 0, 
A„ = Pi(-n(l + a)), n > 0, 
1 7«+i 



(6.4) 



/3n 



n > 1, 



1 + a 7„ 

and the sequence of polynomials qa = 1, and for n> 1 
(6.5) g„(x)=C'^(x)+/3„fc:L^i(a;), 

where is the Krawtchouk polynomial i6.1\) (a ^ 0, — Ij. Consider the second 

order difference operator Da,N 116. 2\) with respect to which the Krawtchouk polyno- 
mials are eigenf unctions. Write finally D for the difference operator of order 2k + 2 
and genre {—k — 1, fc + 1) 

1 



D^Pl{Da,N) + 



1 



-yP2{Da,N)- 



Then D{qn) = A„g„, 71 > 0. 



For a convenient choice of the polynomial Pi, the polynomials (g„)„ in the previ- 
ous Lemma turn out to be orthogonal with respect to a moment functional. Indeed, 
for k a positive integer, a 7^ 0, —1 and N =^ 0, —1, • • • , —k let pk,a,N be the measure 
defined by 

(6.6) Pk.a.N ^ ix+l)---{x + k)pa,N+k+l{x -f fc + 1) 

where Pa.N is the orthogonalizing moment functional for the Krawtchouk polyno- 
mials {k'^-'^)n. For N a positive integer, we have that 



(6.7) Pk,a,N 



(1 



]N+k 



{-ifk\s^k-i + J2 



T{N + k + l)a 



(x + fc + l)r(iV 



x+k+l 



To simplify the notation we remove the dependence of a, N and k and write p = 

Pk,a,N- 

The proofs of the following theorems are similar to the ones in the previous 
Sections and are omitted. 
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Imkr 



eigkrk 



def delkr 



defbetkr 



def qkr 



defplkr 



defp2kr 



ImkrS 



Lemma 6.3. For fc > 1, a ^ 0, -1, iV e 



and n>Q, we have 

a,-N , 



(1 + a)" 



Theorem 6.4. For k > 1, let a and N be real numbers satisfying a 7^ 0, —1 and 



kfj ^ {—n) 0, n > 1, where fc^' ^ , k > 1, are Krawtchouk polynomials (see 



-k 



n 7„+i 



I16.1\) ). We define the sequences of numbers (A„)„, (7n)n>i and (/3„)n>i by 

(6.8) A„ - '^-^^-^+1^ 

(6.9) 7„ 

(6.10) = 

1 + a 7„ 

and the sequence of polynomials {qn)n by qo = 1, and 

(6.11) qn{x) = C"^!^) + /3nC-^i(a;), n>l. 

If N is not a positive integer, then the polynomials {qn)n o-re orthogonal with respect 
to the moment functional p 116.6]) . If N is a positive integer, then the polynomials 
{qn)o<n<N-i o,i"e orthogonal with respect to the measure p \6.'T^ . Moreover, consider 
the difference operator of order 2fc + 2 and genre (— fc — 1, fc + 1) defined by 

D - Pl{Da.N) + -^^P2{Da,N), 

1 + a 

where Da,N is the second order difference operator for the Krawtchouk polynomials 
(/c°'^)„ (see \6.2(l . and Pi and P2 are the polynomials of degrees k + 1 and k, 
respectively, defined by 



(6.12) 
(6.13) 

Then D(g„) = A„g„. 



Pi(x) 



P2ix) = kl 



a-N+1 
k+1 

-N 



Conjecture 2 in the Introduction for the Krawtchouk weight pa,N and the finite 
set F = {1, 2, • • • , fc} can be deduced as a Corollary of the previous Theorem: 

Corollary 6.5. For a positive integer k, write F = {1,2,- •• , fc} and ^f{x) — 
Wj^F^^ ~ j)- Assume that a ^ 0,-1 and ~ ~^{—n) ^ 0, n > 1. Then, 

the orthogonal polynomials with respect to the moment functional pFPa,N i6.6\) are 
common eigenfunctions of a {2k + 2)-order difference operator of genre (— fc — l,fc + 
1), where Pa.N 'is the moment functional for the Krawtchouk polynomials. 



For a + c- N ^ -1,-2,- 
polynomials defined by 



defhap I (7.1) 



'a. C.N 



(x) 



E 



7. HAHN CASE 

we write (ft."'^'^)Ti for the sequence of monic Hahn 

{-n)j{-x)j{l -N+j)n-j{c + j)n-j 

{n + a + c-N + j)n-jjl 



(we have taken a slightly different normalization from the one used in |35| . pp. 
30-53, from where the next formulas can be easily derived; for more details see [25] . 
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deflia3F2 



pp, 204-211). Hahn polynomials can be defined using the hypergeometric function 
3^2: 

,a.c,N( N ^ {c)nil-N)n 'A + a + c- N)j 

" ^' {n + a + c-N)^f^^ {1 -NUc), jl 

(c)„(l - N)„ ( -n, -X, n + a + c- N 



sodeha 



trha 



^^•^^ - {n + a + c-N)^'^\ ' c,l-A^ 

Hahn polynomials are eigenfunctions of the following second order difference oper- 
ator 

(7.3) Da,c,N = x{x - a)5-i + {x + c){x - N + l)Si 

+ [-2x^ + {a-c + N -l)x + a + N{c - 1) - l]5o, 

Da.,c,N{K'''^) = (n + l)(n + a + c-N- l)/i^'"^^, n > 0. 

They satisfy the following three term recurrence formula {h"_'i''^ — 0) 

(7.4) xh„ = hn+i + Khn -I- c„/i„_i, n > 
where 

^ _ c(7V - 1) (a -I- c - TV - 1) -I- 71(0; - c + N - l){n + a + c- N) 
" ~ {2n + a + c- N - l)(2n -ha-Hc-7V + l) ' 

_ n(iV - n)(n + a + c - TV - l)(n + a - iV)(n + c - l)(?i + a + c - 1) 
~ (2n -I- a 4- c - TV - 2)(2n + a + c - AT _ l)2(2n + a + c - A^) 

(to simplify the notation we remove the parameters in some formulas). 

Assume that a-|-c — 7V-|-l,a — A^-fl,a-|-c, ct^O, —1, —2, • • • . If, in addition, 
N is not a positive integer, then the Hahn polynomials are always orthogonal with 
respect to a moment functional Pa,c,N, which we normalize by taking 

, r(a + 1 - jV)r(c)r(a -f e) 

^> - Tia + c+l-N) ■ 

When A'^ is a positive integer, the first N Hahn polynomials are orthogonal with 
respect to the Hahn measure 

(7-5) P.,c,N = T{N) 

(which it is positive when a > A'^ — 1 and c > 0) . 
We also need the so-called (monic) dual Hahn polynomials. 

k 

defdualhap (7.6) /i^' (x) = '-Sj^N-a-dx), 

where Sj^N-a-c, j > 0, are the polynomials defined by (|2.6p . that is, so^N-a-c = 1 
and for j > 1 

i-i 

Sj^N-a^cix) = {-ly Y]_\-^ + ~ ^ - C - «)]■ 

i=0 



haw 
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defdha3F2 



osigj 



ulll 



1721 



1722 



Using Lemma 12.21 we can also write the dual Hahn polynomials in terms of the 
hypergeometric function 3F2 



(7.7) 



{x{x + a + c-N)) = (c)fe(l - N)k J2 
= ic)k{l - N)k3F2 



il-NUc)jj\ 

-k, —X, X + a + c — N 



(|7.2p and (|7.7I) give the duality of the Hahn and dual Hahn polynomials with respect 
to the sequences n and n{n + a + c — N), that is, for k,n > 

(c)„(l - 7V)„/i*'"'^'^(n(n + a + c-N))^ (c)fe(l -N)k{n + a + c- N),,hZ'-'^ {k). 

The computations in the Hahn case are technically more involved because the 
eigenvalue sequence (0„)„ (|7.3p for the Hahn polynomials (with respect to its second 
order difference operator) is quadratic in n. We first need some notation. 

We write (cr„)„ and (6'„)„ for the sequences 



(7.8) 
osigj2 (7.9) 



0n = in + l){n + a + c~ N 
cr„ = 2n + a + c - iV - 2, 



((cn)n is one of sequences which define the first and four P-operators in Lemma 
7.21 below). For j > 0, we also consider the sequence 



Uj (n) = {n)j{—n — a — c + N + 2)j, 
and the polynomials tq = 1 and for j > 1 



defrj (7.10) 



rjix) = {-ly Yl[x + i{a + c- N -2-i)]. 



i=0 



Since the polynomial rj has degree just j, all polynomial P G P can be written as 

-f (2;) = X^^lo'^'' '^j'^'ji^) fo'" certain numbers Wj, j = 0, • • • ,deg(P). 
We will need the following technical Lemma. 

Lemma 7.1. Let a, c, A'^ G M be real numbers. Then for j,n > 



(7.11) 

(7.12) anUj{n) + an+iUj{n+ 1) = -2 



Uj+i{n + 1) - Uj+i{n) 
J + 1 



+ {-a -c + N + 2{j + l))(uj(n + 1) - Uj{n)). 

Proof. (|7.1ip is just Lemma [2.21 for x = —n and u — a + c~N — 2. 
(|7.12p of the Lemma can be checked by a direct computation. 



□ 



Proceeding as in the previous Sections, we have identified four P-operators of 
type 2 for Hahn polynomials (the proof is omitted) . 

ITha Lemma 7.2. For a, c, N satisfying a + c — N ^ —1, —2, • • • , consider the Hahn 
polynomials (/i^''^'^)n 17. ij). Then, the operators T>i, i = 1,2,3,4 defined by 113. 9\) 
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th5 . lha 



from the sequences (n>0) 

n{N - n){n + a - N) 



En,! 
Era, 2 
En, 3 
En, 4 



{2n + a + c- N - l)(2n + a + c- N ~2)' 
n{n + a — N){n + a + c — 1) 

{2n + a + c- N - \){2n + a + c- iV-2)' 
-?i(7V - + c - 1) 

{2n + a + c- N - \){2n + a + c - N ~ 2)' 
—n{n + c — l)(n + a + c — 1) 



a„,i ^2n + a + c~ N -2, 
o-n,2 = -(27i + a + c - iV - 2), 
cr„,3 = -{2n + a + c- N -2), 
(Tn 4 = 2n + a + c - N - 2, 



{2n + a + c- N - l)(2n + a + c- iV-2)' 

are V-operators for {h'^''^'^)n and the algebra A fl.lO\) of difference operators with 
polynomial coefficients. More precisely 

a + c — N 
V, = {-X + TV - 1)A - ^—^ /, 

^ / N„ a + c — N ^ 
V2 = {x- a)V + /, 

^ „ a + c- N ^ 
P3 - a;V + /, 

„ / N . a + c — N 
Vi = -(.T + c)A ^ /. 

Each one of the I?-operators displayed in the previous Lemma together with 
Lemma 13.11 gives rise to the corresponding class of polynomials satisfying higher 
order difference equations. We only proceed with the first and second 2?-operators 
because due to the symmetries of Hahn polynomials, the examples generated from 
the first and third operators are essentially the same, as well as the examples 
generated from the second and the fourth operators. 

Theorem 7.3. Let a, c, N be real numbers satisfying a + c — N ^ —1, —2, • • • . Let 

P2 be an arbitrary polynomial of degree k > 1 which we write in the form 

k 

P2{x) ^^Wjrj{x), 

for certain numbers Wj, j = 0, • • • Consider also the polynomial Pi of degree 
k + 1 defined by 

k 



Pi{x) = {a + c-N)P2{x) + 2{x-a~c + N 



w 



J + 1 



Consider next de sequences of numbers (7„)„, 
by 

k 

= ^WjU.j{n), n>0, 



i)„, « = 1,2, and (A„)„ defined 



obnj (7.13) 



obnj2 (7.14) 



/3„,i = 

/3n,2 = 



i{N ~ n){n + a- N) 



In+l 



{2n + a + c- N -l){2n + a + c- N -2.) 7„ ' 
n[n + a - N){n + a + c - I) In+i 
{2n + a + c~ N -l){2n + a + c- N -2.) 7„ ' 



n > 1, 
n > 1, 
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inj (7.15) A„ = , n > 1, Ao = 



qiinoj 



alPl 



2 2 

where (On)n and ((t„)„ are defined by \7.8^ and |y.g[ j, and we implicitly assume 
that 7n 7^ 0, n > 1. Define finally the sequences of polynomials qo^i = 1, and for 
n>l 

(7.16) g„,, = /iJJ^^'^ + /3n,^h'^'^f, * = 1, 2, 

where h"''^'^ is the Hahn polynomial {7.1^ . Consider the second order difference 
operator Da^c.N 1^7. 3\ ) with respect to which the polynomials (ft-^'^'^)n are eigen- 
functions. Write finally Di, i = 1,2, for the difference operators of order 2k + 2 
and genre {—k — 1, fc + 1) 

Di = -Pi{Dc,^c.n) + f ^ + {~X + N- 1)A j P2(i?a,c,iv), 

1 f a + c- N \ 

D2 = -PliD^^cM) - ( ^ + (X - a)Vj P2(i?ax,Ar). 

Then Di{qnA) = Kqn,i, n > 0, i = 1, 2. 

Proof. The theorem is a straightforward consequence of the Lemma [3.21 We have 
just to identify who the main characters are in this example. 
Write 

_ n{N - n){n + a - N) _ 

^ {2n + a + c-N- l)(2n + a + c - iV - 2) ' '^".1 = '^"' 
n(n + a — N){n + a + c — 1) 
(2n + a + c-iV-l)(2n + a + c-iV-2)' ^".2 = 

Lemma 17.21 gives that the sequences (£„.i)„ and (tT„_i)„, i = 1,2, define two 2?- 
operators of type 2 p.9p for the polynomials {h"''^'^)n and that 

2?i = - 2 ^ +i-x + N~ 1)A, 

_ a + c- iV , 

V2 = / +ix- a)V. 

Consider now the polynomial 

(7.17) Q(x) = f-±-r,+,{x) + {-a-c + N + 2(j + l))r,(x)' 

A straightforward computation using (j7.10p gives that Q — Pi. 
According to the Lemma we have just to check that 

(1) Jn+l=P2{0n),n>0. 

(2) A„ - Xn-i = cr„7„, n>l. 

(3) A„ + A„_i = Pi{On-i), n > 1, where 0„ are the eigenvalues of {h"''^'^)n 
with respect to Da,c,N- 

Formula (1) is a direct consequence of (|7.1ip in Lemma mi 
Using the definition of A„, one can easily see that for n > 2 conditions 2 and 3 
above are equivalent to 

-Pl(6'«) - Pl{dn~l) _ CFnln + Cr„+l7„+l 
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dhp 



sodedh 



f odedh 



But this is a direct consequence of (|7.12p in Lemma 17.11 and the formula (|7.17p for 
Pi- 

Conditions 2 and 3 above for n = 1 can be checked by a straightforward compu- 
tation. □ 

7.1. Hahn I. For a convenient choice of the polynomial P2, the polynomials {qn,i)n 
in the Theorem 17 . 31 turn out to be orthogonal with respect to a moment functional. 
To identify that polynomial P2, we need the family of dual Hahn polynomials 



*,N+c-l,2-c,a+c-l 



n.fe 



(7.18) hl,{x) 
hl^{x{x + 2 + N 



E 

3=0 



Using (fTHI) and (fTTT]) we get 

(-fc),(2 -a-c + j)k-ji2 - c + j)k-, 



rj{x), 



c)) 



= {2 - a - c)k{2 - c)k3F2 



-k, -x,x + 2 + N - a- 
2-a-c,2-c 



1 



This coincides with the hypergeometric representation for the (monic) dual Hahn 
polynomials given in [23], p. 208 (with the notation of [23] . 



hlJX{x)) = (2 



c)fe(2 - c)kRkiXix), 1 - a - c, TV, c - 2) 



where A(x) = x{x + 2 + N — a — c)). They satisfy the following second order 
difference equation in the lattice x{x + 2 + N — a — c) (see pHj pp. 209) 

(7.19) r(.x)S-i(/i^fe) - ir{x) + .s(x))So(/it,fe) + s{x)Siihlk) = /cu(x)So(/it,fe), 

where Si{p) ^ p{{x + l){x + 1 + 2 + N - a - c)), I G Z, and 

r{x) = -x{x + N){x -a + N){2x - a- c + N + 3), 

(7.20) 

s{x) = {x-a-c + N + 2){x -a-c + 2){x - c + 2)(2a: - a - c + N +\), 
u{x) = (2x - a - c + N + l){2x - a - c + N + 2){2x ~ a - c + N + i). 

They also satisfy the following first order difference equation ([23] pp. 210): if we 
write A = Si — So, then 



(7.21) K{hl'^^)^k{2x~a-c + N + i)h*{^-_\{x{x + i + N ■ 

For k a positive integer, and a, c, G M satisfying 



c)). 



conacN 



thaw 



thaw2 



(7.22) a + c-7V + l,Q:-7V+l,a + c-A:-l,c-fc-ly^ 0,-1, -2, 
let Pk,a,c.N be the moment functional defined by 

(7.23) Pk,a,c,N = (x + c - 1) • ■ • (a; + c - k)pa,c~k~i,N 

where Pa,c,N is the orthogonalizing moment functional for the Hahn polynomials 
(/i"''^'^)„. When iV is a positive integer we have that 

N-l 



(7.24) 



Pk,a,c,N 



= r(iv)E 



r(a — x)T{x + c) 



-1 (x + c - 

2:— ^ 



i)r(7v- 



X }x\ 



(only the first N Hahn polynomials are then orthogonal) . To simplify the notation 
we remove the dependence of a, c, N and k and write p = pk,a,c,N- 
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In order to find the orthogonal polynomials with respect to pk,a,c,N we need the 
following Lemma. 

Lemma 7.4. For k a positive integer, let a, c, N be real numbers satisfying {7.22^ . 
Then, we have 

(-l)"n!(iV - n)„r(c - fc - l)r(a + n-N+ l)V{a + c-k- l)hl,{e,,) 
{Pk,^,c,N, K" )- r{a + c-N + 2n) 

where {9n) is the sequence defined by |7.(g[ j. 

Proof. We proceed in an analogous way to the Lemmas 14 . 2 1 or [531 For n > write 

_ (-l)"n!(7V - n)„r(c - fc - l)r(a + n-N+ l)T{a + c-k~ l)/it_fe(6l„) 
^" ~ T{a + c- N + 2n) ' 

The three term recurrence relation ()7.4p for {h'^''^'^ )n and the normalization for the 
Hahn weight Pa,c,N give that 

I 1 .X , „ , X. r(c)r(Q; + 1 - Ar)r(Q; + c) 

Ultdv I (7.26) ^n+l + {bn + C-k-l)^n+CnCn-l=0, n>l, 

where 6„ and c„ are the recurrence coefhcients of the Hahn polynomials (/i"''^'^)n. 
We now prove that also 

I 1 > /, , r(c)r(a + 1 - 7V)r(a + c) „ 

mecmeSa 7.27 Cl + fa + c - fc - 1 Co - ^ ') ^ ^' ^ = 0, 

mecme3 | (7.28) C«+l + (&« + c - fc - 1)C„ + c„C«-l = 0, n > 1. 

After straightforward computations taking into account the definition of (n , (|7.28p 
is equivalent to 

r(-n - l)/ij,fe(0„+i)-(r(-n - 1) + s{-n - 1) + ku{~n - l))hlM 

+ si-n - l)hlk{en-i) ^ 

where the polynomials r, s and u are given by (j7.20p . But this follows just by 
writing x — —n — 1 in the second order difference equation (|7.19p for the dual Hahn 
polynomials (ft.*^fe)fe. (|7.27p follows in a similar way. 

Since the sequences {^n)n and (Cn)n satisfy the same recurrence relation, it is 
enough to prove that = Co- If we write 

Vk.c = Co = {Pk,a,c.N, 1) = {Pa.c.N, [x + C - I) ■ ■ ■ [x + C - fc)), 

and proceed as in Lemma 15.31 we find that 

r(c - l)r(a - AT + i)r(a + c - 1) 

rik^c - fc7?fe-l,c-l = ^:r^ ■ TT. • 

On the other hand, if we write 

r(c - fc - l)r(a - TV + l)r(a + c - fc - l)h*{l{a + c-N -1) 



Tk,c = Co 



T{a + c-N) 
a simple computation using (|7.2ip gives 

_ r(c - i)r(a - TV + i)r(a + c - 1) 

I [a + c — N) 
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lm61ha 



defplha 
defp2ha 



eighak 



def delha 



defbetha 



def qha 



cha555 



From where we easily get = Co 



□ 



We are now ready to prove that the orthogonal polynomials with respect to 
Pk,a,c,N (|7.23p are a particular case of Theorem 17.31 



Theorem 7.5. For k > I, let a, c, TV be real numbers satisfying and /ij f^{n{n- 

2 — iV + a + c)) ^ 0, n > 0, where h\ j^, k > 1, are the dual Hahn polynomials ^7.18^ . 
We then define the polynomials P2 and Pi, of degrees k and fc + 1, respectively, by 

(7.29) P2{x)^hl^{x), 

(7.30) Pi{x) = [a + c- N)hlj,{x) + 2{x^a-c + N +1) 



k 

xE 

J=0 



(-fc)j (2 - a - c + j)k-, (2 - c + j)k-, 
We also define the sequences of numbers (A„)„, (7n)n>i o,nd (/3„)„>i by 

CTnln + Pi 



,{x). 



(7.31) 
(7.32) 

(7.33) 



-1) X PM~<JiP2i9o) 

, n>l, Ao = 



In = hl^{n(n -2-N + a + c)), 

n{N - n){n + a - N) 



2 

7«+i 



{2n + a + c- N -l){2n + a + c- N -2) j„ 
and the sequence of polynomials {qn)n by qa = 1, and 
(7.34) (7n(x)=/.r'^(x)+/3„C-f(x), n>l, 



where {9n)n o,nd (cr„)„ are defined by \7.^ and \7.9^ . If N is not a positive integer, 
then the polynomials (g„)„ are orthogonal with respect to the moment functional p 
\7.2S^ . If N is a positive integer, then the polynomials {qn)o<n<N o,i"£ orthogo- 
nal with respect to the measure p {7.24^ . Moreover, if we consider the difference 
operator of order 2fc + 2 and genre {—k — l,fc + 1) defined by 



D = -Pi(Da,,,,N) 



N 



+ i-x + N -l)A]P2iD^.c^N) 



where D^^cN "is the second order difference operator for the Hahn polynomials |7.i5[ ), 
then D{qn) ^ A„g„. 

Proof. The first part of the Theorem can be proved in an analogous way to Theorem 

1131 

The second part is the particular case of Theorem 17.51 for z = 1 and 

(-fc)j (2 - a - c + (2 - c + j)k-j 



□ 



Conjecture 1 in the Introduction for the Hahn weight Pa,c,N and the finite set 
F = {1, 2, • • • ,k} can be deduced as a Corollary of the previous Theorem: 

Corollary 7.6. For a positive integer k, write F = {1, 2, • • ■ , fc} and p^c-pix) = 
Yijepi^ + c + j). Let a, c, N be real numbers satisfying that a + c — N-\-k + 2,a — 
N + l,a + c,c ^ 0, —I, —2, ■ ■ ■ and that 

hlf+^'+^'^'inin + a + c + k-N-l))^0. 
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Then, the orthogonal polynomials with respect to the measure p-c-FPa,c.N are com- 
mon eigenf unctions of a {2k + 2)-order difference operator of genre (— fc — 1, fc + 1). 

7.2. Hahn II. There is other choice of the polynomial P2 for which the polynomials 
{111,2)11 in the Theorem [731 are also orthogonal with respect to a moment functional. 
Indeed, write {h^ f.)k for the family of dual Hahn polynomials hl^ j. — /i^'""'^""^'"^. 
Using ILH) and (jfjl) we get 



dhp2 (7.35) 
(7.36) 



{-k)j{2-c + j)k-AN + l+])k-, 



rj{x), 



h*2^{x{x + 2 + N ~ a ~ c)) 



,^ \ ,i.T -,\ / -k, -x,x + 2 + N - a — c ^ 
= {2~cUN + l),,F2[ + i ;l 



This coincides with the hypergeometric representation for the dual Hahn polyno- 
mials given in 23 , p. 208 (with the notation of |23) . 

KAHx)) = (2 - c)k{N + l)kRk{X{x), l~c,N~a,-N-l) 

where X{x) ~ x{x + 2 + N ~ a — c)). 



conacN2 



thawb 



thawb2 



Im61ha2 



For k a positive integer, and a,c, N ^ M. satisfying 

(7.37) a + c- N +l,a - N +1, a + c,c - k -1 ^0,-1,-2, - ■■ , 
let pk,a,c.N be the measure defined by 

(7.38) Pk,a,c,N = {x+l)---{x + k)pa+k+l,c-k-l,N+k+l{x + k+ 1). 

When is a positive integer we have 



(7.39) 



Pfc,a,c,w = {-lfk\T{a + + l)r(c - fc - \)5-k^i 

^ ^ T{a — x)T{x + c) 



r(iv + fc 



x=0 



(a; + fc + l)r(7V - a;)a;! 



To simplify the notation we remove the dependence of a, c, N and fc and write 

P = Pk,a,c,N- 

The proofs of the following Theorems are similar to that of the previous Sections 
and are omitted. 



Lemma 7.7. For fc a positive integer, let a,c, N be real numbers satisfying {7.31^ . 
Then, we have 



{Pk,a,c,Nih"'''' ) 



(-l)"+'=n!r(c - fc - l)T{a + c + n)T{a + l + n- N)h*j^{en) 



r{a + c-N + 2n) 
where (0„) is the sequence defined by \7.8^ . 



Theorem 7.8. For k > I, let a, c, N be real numbers satisfying |y.,?7| ) and /ij fe('T^(?^ 
2 — N + a + c)) 0, n > 0, where h^ k > I, are the dual Hahn polynomials Iji7.35\ ). 
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defp21ia2 



eigmek2 



defdelha2 



defbetha2 



hah aha 



We then define the polynomials P2 and Pi, of degrees k and fc + 1, respectively, by 

(7.40) P2ix) = hl.ix), 

(7.41) Pi{x) = (a + c - N)hl,,{x) + 2{x - a - c + N + I) 

^ (-fc),(2 -c + j)k^,{N+l+j)k^, 



We also define the sequences of numbers (A„)„, ijn)n>i and (/3n)n>i by 



(7.42) 
(7.43) 

(7.44) 



A. 



, n>l, Ao = ; 



7n KA'^i'^ -2- N + a + c)), 

n{n + a — N){n + a + c — 1) 

Pn — 



ln+1 



{2n + a + c- N -l){2n + a + c- N -2) 7„ ' 
and the sequence of polynomials {qn)n by qo — 1, and 



defqha2 (7.45) 



qn{x)^hr'^{x)+P„K^lf{x), n>l, 



N , 



where {9n)n and (cr„)„ are defined by Ii7.8\ l and |y.ff[ j. If N is not a positive integer, 
then the polynomials {qn)n are orthogonal with respect to the moment functional p 
^7.38^ . If N is a positive integer, then the polynomials {qn)a<n<N are orthogo- 
nal with respect to the measure p \7.39^ . Moreover, if we consider the difference 
operator of order 2k + 2 and genre {—k — l,fc + 1) defined by 

D = -Pi{D^,c,n) - i ^ + {x- a)Vj P2{D^^c,n), 

where Da c.N is the second order difference operator for the Hahn polynomials 117. 3\ l, 
then D{q„) = A„(7„. 

Conjecture 2 in the Introduction for the Hahn weight Pa,c,N and the finite set 
F = {1, 2, • • • , fc} can be deduced as a Corollary of the previous Theorem: 

Corollary 7.9. For a positive integer k, write F = {1,2,- •• , fc} and ppix) = 
Yijepi^"^)- ^ be real numbers satisfying a+c—N+k+2, a—N+1, a+c, c ^ 

0, — 1, • • • , and that 



h. 



*,a-k-l,c+k+l,N~k~l 
2,k 



(n{n + a + c- N + k-l))^Q. 



Then, the orthogonal polynomials with respect to the moment functional p FPa,c.N 
are common eigenf unctions of a {2k + 2)-order difference operator of genre (— fc — 
l,fc+l). 



8. Appendix 

In this Appendix, we apply our method to the Laguerre and Jacobi polynomials 
to generate families of orthogonal polynomials which are eigenfunctions of higher 
order differential operators. 
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deflap (8.1) 



Laguerrew 








sodolag 








ITlag 



th5. Hag 



becc 



qimolag 



8.1. Laguerre polynomials. For a G R, we use the standard definition of the 
Laguerre polynomials {L^)n 

{—xy / n + 



J! 



(that and the next formulas can be found in [TT], pp. 188-192). 

When a ^ — 1,— 2,---, they are orthogonal with respect to a measure /i^ — 
liaix)dx^ which it is positive only when a > —1, and then 

(8.2) ^„(x) = a;"e"^, < x. 

Laguerre polynomials are eigenfunctions of the following second order differential 
operator 

2 , 



(8.3) 



Da = X 



1 



^dx J dx 
We first identify a D-operator for the Laguerre polynomials. 



~riL^, n > 0. 



Lemma 8.1. For a G M, let Pn — L", n > 0, the Laguerre polynomials given by 
iS.l]) . Then the operator V defined by h3. ^|) from the sequence £« = —1, n>Q, is a 
V-operator for the Laguerre polynomials and the algebra A 11.11]) . More precisely 
V = d/dx. 

Proof. Since p„ = L", using the definition above of the operator P and the well- 
known formulas (£")' = -L"tl, L^+^ = J2'j=o^n-]' ^^^^ 



n n — 1 



3=0 



dx 



This gives V = d/dx. 



□ 



The previous Lemma together with Lemma 13.11 gives rise to the corresponding 
families of polynomials satisfying higher order differential equations. 

Theorem 8.2. Let Pi be a polynomial of degree k + I, k > 1, and write P2{x) — 
Pi{x~l)—Pi{x) (so that P2 is a polynomial of degree k). We assume that P2{—n) ^ 
0, n > 0, and define the sequences of numbers 

7„+i = P2i~n), n > 0, 
A„-Pi(-n), n>0, 
7«+i 



(8.4) 



I3n = — 



n > 1, 



In 



and the sequence of polynomials qo = 1, and for n>l 
(8.5) g„-L^ + /3„L^_i, 

where L" is the n-th Laguerre polynomial. Write finally D for the differential 
operator of order 2k + 2 



D = Pi{Da 



dx 
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darlag 



14.1 



def dellagnn 



def betlagnn 



def qnlagnn 



kkk 



eiglagk 



where Da is the second order differential operator with respect to which the Laguerre 
polynomials {L")n are eigenfunctions. 
Then D{qn) = Xnqn, n>0. 

The case fc = 1 in the previous Theorem is included in |13) . and for general k is 
implicit in [19j. 

We next consider some important particular cases of the previous theorem. 
For a ^ 0, —1, —2, • • • and M ^ —T{n)/T{n + a), consider the measure pa^M 
defined by 

(8.6) pa^M = c(T^{a)M6a + Ha-iix)dx, < x, 

where Pa{x) is the orthogonalizing function for the Laguerre polynomials (L^)„. 
When a > and A/ > 0, the measure Pa,M is the Krall-Laguerre-Koornwinder 
measure 



aV'^{a)M5o + x"^^e- 



Q<x. 



To simplify the notation we remove the dependence of a and M and write p — Pa,M- 
Using Lemma 12.11 we find explicitly a sequence of orthogonal polynomials with 
respect to p in terms of p„ = L", n > 0. 

Lemma 8.3. Let (7ri)n>i be the sequence of numbers defined by 

T{n + a) 



(8.7) 

Write 



7„ = 1 + M- 



7n+l 



r(n) 



n> 1. 



Then the polynomials defined hy qo = 1 and 

(8.9) qn = LZ+PnLl_^, n>l, 

are orthogonal with respect to p h8. Oj) . 

Proof. In the notation of Lemma 12.11 have \ — v — pa~i{x)dx and /x = 
Pa(x)dx. Hence the formula — J2^=o -^n-j gives a„ — T{a). Since we also 

have p„(0) = ("^"): an straightforward computation gives (|8.7p and (18.81) from 



(|23D . 

The formula 



□ 



is the particular case fc = 1 of formula (3.16) in 



When a is a positive integer, the sequence (7n)n (|8.7p is a polynomial in n, 
and then Theorem 18.21 gives the (2a + 2)-th order differential equation that the 
Krall-Laguerre-Koornwinder polynomials (|8.9p satisfy. 



Theorem 8.4. Let (g„)„ be the sequence of orthogonal polynomials k8. 9fl with re- 
spect to the measure p 18. 6\) . Assume that a is a positive integer. Define the 
sequence of numbers A„ by 

M 

(8.10) A„=n+ -n{n + l)---{n + a), 

a+l 

and the polynomials Pi and P2, of degrees a + l and a, respectively, by 

M 



Piix) 



1 



{-x){-x + l)---{-x + a), 



P^{x) = I + M{~x + l){-x + 2) ■ ■ ■ {-X + a). 
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Consider the differential operator of order 2a + 2 defined by 



pqdlc2 



pqdlcS 



pqdlc 



D 



where Da is the second order differential operator for the Laguerre polynomials 
defined by \8. 3\) . Then 

D{qn) = A„g„. 

Proof. It is enough to check that P2{x) = Pi{x — 1) — Pi{x) and then use Theorem 



□ 

Except for a sign, the eigenvalues (A„)„ (|8.10|) coincide with that given by J. and 
R. Koekoek in [20] ■ Hence, we can conclude that the (2a + 2)-th order differential 
operator —D, where D is displayed in the previous theorem, coincides with the dif- 
ferential operator found by J. and R. Koekoek for the Krall-Laguerre-Koornwinder 
polynomials. This provides a nice expression for Koekoek differential operator as a 
polynomial combination of the Laguerre second order differential operator. 



When a is not equal to fc, the polynomials (q„)„ (j8.5|) in Theorem [521 seem not to 
be orthogonal with respect to a measure. Anyway they enjoy certain orthogonality 
property. Indeed, when -P2(l) 7^ 0, since P2 is a polynomial of degree k, we can 
choose numbers Wj, j = 1, ■ • • , /c, such that 

(8.11) P2i-x)^P2{l) 




On the other hand, if ^2(1) = 0, we can choose numbers Wj, j = 0, 
that 



, fc, such 



.12) 



P2{-X) = 1 



k 
J=0 



x + j 
j 



(in particular wq = —1). 
We then define the polynomial Q by 



.13) 



Q(x) = 




j)jW.jX^ 
-i)jWjX^~3, 



if ^2(1) ^ 0, 
if F2(l) = 0. 



Notice that the polynomial Q has degree at most /c — 1 for P2(l) 7^ 0, but k for 
^2(1) = 0. 

Lemma 8.5. Let P2 be an arbitrary polynomial of degree k satisfying that P2 (~'t-) 7^ 
0, n > 0. For a G M \ {fc, fc — 1, • • • }, consider the polynomials {qn)n defined by 
h8.5\) . If we write 



iimlag (8.14) 



f{x)g{x)^ia-i{x)dx + g{0) / f{x)Q{x)^ia-k-iix)dx, 



where fia is the orthogonalizing measure for the Laguerre polynomials (i^)n and 1 
is the polynomial defined by I18.13\) . then 

(1) {Qn,qj) 0, j = O,--- ,n-l, 

(2) {qn,qn) 7^ 0, and 

(3) {x'^+'^qn, q-j) ^0, n>k + 2 andj^O,--- , n - fc - 2. 
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Proof. It is enough to prove that (q„, x-') = 0, j = 0, • • • , n — 1, (g„, a;") ^ 0, and 
{x''+'^q„,x^) = 0, 71 > A; + 2 and j = 0, • • • ,n-k-2. 

Since qn is a hnear combination of L" and L^_i, we automaticahy have that 

{qn,x^) ==0, j = !,■■■ ,n-l, 

{x^+^qn,x^) =0, n>k + 2 and j = 0, • ■ • ,n - fc - 2, 

whatever the polynomial Q and the numbers /3n, n > 0, are. In the same way, we 
have for n > 1 

= y qn(x)x"-^fi^{x)dx = {-ir-\n-l)l/3jL^^_,\\2^Q. 

For j = 0, we proceed as follows. Using the expansion — X]j=o ^^~jf^^n-j 
(see [ri], p. 192, (39)) we get for ^ > 1 

L>„_i {x)dx = r{a~l + l) Y_ ^ 

We first prove the Lemma when ^2(1) ^ 0. Using the definition of Q (|8.13p and 
taking into account that x-'/Iq = Ha+j, we get 

n + j\ f n + j — 1 



w 



3 J \ 3 



(g„,l) =r(a) |^(l + /3«) + E 
Equating ((?„, 1) = 0, we get (using (|8.1ip ) 

1 + E,ti-.(T) _ P^{-n 



= — 



l + E.ti^.CT') ^2(-(n-l))' 
which it is precisely the expression for /?„ in Theorem 18.21 
Finally, 

(i,i) = r(a) (^i + J2w}j =r(a)P2(o)/P2(i)^o. 

When ^2(1) — 0, the proof is similar. □ 

The orthogonality properties above for the polynomials {qn)n imply that they 
have to satisfy a (2k + 3)-term recurrence relation. 

Corollary 8.6. Assume that a 7^ fc, fc — 1, • • • , then the polynomials {qn)n defined 
in Theorem \8.2\ satisfy a (2k + 3) -term recurrence relation of the form 

k+l 

x''^^qn(x) = ^ Cn.,jqn+j(x). 
]=-k-l 

Proof. Write x'^'^^qn ~ J2j=-n'^n,jqn+j- Then, part (3) of the previous Lemma 
gives 

0= {x''+'^qn,qo) ^ Cn-n{qo,qo)- 

Using part (2) of the previous Lemma we get c„^_„ = 0. We now proceed by 
induction on m. Assume c„._„+j = 0, j = 0, • • • , m, m < n — fc — 2. Proceeding as 
before, we find 

~ {x''^^ qn, q-n+m+l) = Cn,-n+m+l{q-n+m+l, q-n+m+l) , 
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def jac 



jacw 



sodo j ac 



IT jac 



from where we deduce that also o 



= 0. 



8.2. Jacobi polynomials. For a, /3 G M, a + /? 7^ — 1, —2, ■ 
definition of the Jacobi polynomials {Pn'^)n 



□ 

we use the standard 



.15) 



n + a 
j 



n + 
n- j 



j=o 

(that and the next formulas can be found in [TT], pp. 169-173). 

For a, /3, a + /3 ^ —1, —2, ■ • • , they are orthogonal with respect to a measure 
l-J'a.p — fJ'a,i3ix)dx, which it is positive only when a, (3 > —1, and then 

(8.16) ^ia.,pix)^{l~x)°'{l + xf, -l<x<l. 

They are eigenfunctions of the following second order differential operator 

2 



.17) 



+(/3-a-(a + /3 + 2).)A, 



dx 



dx 



Da^APn") = + a + l3 + 1)P,^^^ n > 0. 



We omit the proofs in this Section, since they are essentially the same as in the 
previous Sections. 

We first construct two P-operators for the Jacobi polynomials. 

Lemma 8.7. For a, /? G M satisfying a + (3 ^ —1, —2, • • ■ , let pn = Pn'^ ) n >Q, 
the Jacobi polynomials given by h8.15\) . Then the operators Di and I?2 defined by 
from the sequences (n > 0) 

n + a 



£n.l 



£n.2 



n + a + P 
n + /3 
n + a + 



cr„^i = 2n + a + /? - 1, 



-(2n + a + /3- 1), 



are "D- operator for the Jacobi polynomials and the algebra A 111.11]} . More precisely 



a + f3 + l 



a + B + 1 , , d 

To establish the main results in this Section, we need some notation. 
We write (cr„)„ and (&„)„ for the sequences 

e„ = -n{n + a + P + l), 
an = 2n + a + P - 1. 

We recall that {9n)n is just the eigenvalue sequence for the Jacobi polynomials P^' 
(|8.15p . For j > 0, we also consider the sequence 

Uj{n) {n + a)j{n + (3 - 

and the polynomials rg = 1 and for j > 1 



r-iix) ]^[-a; + (a + i + l)(/3-i)]. 



i=0 
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th5 . Ijac 



Since the polynomial rj has degree just j, all polynomial P G P can be written 

as P{x) = J2'j=o^^ '^j'f'ji^) foi" certain numbers Wj, j = 0, - ■ ■ ,deg(P). 

Lemma 13.21 and the two P-operators for the Jacobi polynomials found in Lemma 
18.71 automatically produce a large class of polynomials satisfying higher order dif- 
ferential equations. We only proceed with the first 2?-operator because due to the 
symmetry of the Jacobi polynomials with respect to the parameters a and /3 and 
the endpoints of its support x = 1 and x = —1 the examples generated by both 
©-operators are essentially the same. 

Theorem 8.8. Let P2 he an arbitrary polynomial of degree k > 1 which we write 
in the form 

k 

P2{x) ^^Wjrj{x), 

for certain numbers Wj, j = 0, • • ■ ,k. Consider also the polynomial Pi of degree 
k + 1 defined by 



3=0 



Piix) = ( 7TT'^^+'^'''^ + {a-l3 + 2j + l)r,{x) 



Consider next the sequences of numbers (7„)„, o,nd (A„)„ defined by 



In ='YwjUj{n), n>0 

3=0 

n + a 7„+i 



/3n = 
An = 



n + a + 13 7„ 

0-„7„ -I- Pl{On-l) 



n > 1, 



n > 1, Ao = 



Pi(0)-(a + /? + l)P2(0) 



qnnoj ac 



where we implicitly assume that Jn ^ 0, n > 0. Define finally the sequence of 
polynomials qo = 1, and for n > 1 



(8.18) 



pa, 13 I o pa, (3 



where P"'^ is the Jacobi polynomial Ii8.15\) (a + (3 ^ —1, —2, ■ ■ ■ ). Write finally D 
for the differential operator of order 2k + 2 



D 



-Pl{Da„ 



a + + 1 



il-x) 



dx 



where D^.p is the second order differential operator with respect to which the poly- 
nomials {Pn'^)n o-TC cigcnfunctions (see 18.1'1\) ). 
Then D{qn) = XnQn, n > 0. 

The case k — 1 in the previous Theorem is included in [13], and for general k is 
implicit in [T8] . 

We next consider some important particular cases of the previous Theorem. 
For a,P,a + l3^ -1, -2, • • • and M ^ -r(n)r(n + a)/{T{n + a + (3)T{n + /3)), 
consider the measure Pa,p,M defined by 



dcirj ac 



.19) 



Pa,l3,M 



_ Oa+0 



/3T^P)M6.i + fi^,p.i{x)dx, 
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def deljac 



def qnj ac 



eigjack 



where iJ.a.p{x) is the orthogonahzmg function for the Jacobi polynomials (P"'^)„. 
When a > — > and AI > 0, the measure Pa,0,M is the Krall-Jacobi- 
Koornwinder measure 



2°+'^/3r2(/3)M(5_i + (1 - + xf-\ 



-1 < a; < 1. 



To simplify the notation we remove the dependence of a, /3 and M and write 

P = Pa,f3,M- 

Using Lemma 12.11 we find explicitly a sequence of orthogonal polynomials with 
respect to p in terms of p„ = P"'^ : n > 0. 

Lemma 8.9. Let (7„)„>i be the sequence of numbers defined by 

T{n + a + p) r{n + l3) 



(8.20) 
Write 



7„ = 1 + M 



T{n + a) 
n + a 7„+i 



n + a + f3 7„ 
1 and 



Tin) 



n> I. 



n> I, 



Then the polynomials go 
(8.21) g„, = P° 

are orthogonal with respect to p i8.19\) . 

The formula ()8.9p is the particular case A: = 1 of formula (3.10) in [18] . 

When /3 is a positive integer, the sequence {jn)n (|8.20p is a polynomial in n, 
and then Theorem 18.21 gives the (2a + 2)-th order differential equation that the 
Krall-Jacobi-Koornwinder polynomials (|8.2ip satisfy. 

Theorem 8.10. Let ((/„)„ be the sequence of polynomials 118. 2 1]) orthogonal with 
respect to the measure p Ii8.19\} . Assume that 13 is a positive integer. Define the 
sequence of numbers A„ by 



..22) 



\n = {n + a + f3)[n + M{n + a)^(n)^ 1 + 



n - 1 



and the polynomials Pi and P2, of degrees /3 + 1 and P, respectively, by 

2x 



Pi{x) = a/3 + (a + l)(/3 + 1) - 2a; + M 



13 + 1 
ip 

+ a + l3 + l 



a/3, 



i=0 



0-1 



P2{x) = l + M\[[-x+{a + i + l)(/3 - i)]. 

i=0 

Consider the differential operator of order 2/3 + 2 defined by 
1 



D^-Pi{D„, 



a+£+l d 



where Da, (3 is the second order differential operator defined by ^8.17^ . Then 

D{qn) = A„g„. 
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Comparing the eigenvalues (A„)„ (|8.22p with that given by Zhedanov in [38j . 
(8.14), we can conclude that the differential operator found by Zhedanov for the 
Krall-Jacobi-Koornwinder polynomials orthogonal with respect to (jl.l6p . = 0, 
coincides (up to normalization) with the (2/3 + 2)-th order differential operator 
displayed in the previous theorem. 

When /3 is not equal to fc, the polynomials (g„)„ (|8.18p in Theorem 18.81 seem not 
to be orthogonal with respect to a measure. Proceeding in the same way as in the 
previous Section (see Lemma l8.5|) one can see that they enjoy similar orthogonality 
properties with respect to the nonsymmetric bilinear form 

{f,g) = J f{x)g{x)fia,i3-i{x)dx+g{-l)J f{x)Q{x)^a.i3-k-i{x)dx, 

where fj,a,fj is the orthogonalizing measure for the Jacobi polynomials {P"'^)n, 
and Q is certain polynomial of degree at most fc — 1 or fc depending on P2 . These 
orthogonality properties for the polynomials {qn)n imply that they satisfy a (2/c+3)- 
term recurrence relation of the form 

k+l 

[1 + X)''^^ qn{x) ^ ^ Cn,jqn+j{x). 

3=-k-l 

References 

[1] R. Alvarez-Nodarse and A.J. Duran, Using "D-operators to construct orthogonal polyno- 
mials satisfying higher order g-difference equations, in preparation. 

[2] H. Bavinck and H. van Haeringen, Difference equations for generalizations of Meixner 
polynomials, J. Math. Anal. Appl. 184 (1994), 453-463. 

[3] H. Bavinck and R. Koekoek, On a difference equation for generalizations of Charlier 
polynomials, J. Approx. Theory 81 (1995), 195-206. 

[4] T. Chihara, An introduction to orthogonal polynomials, Gordon and Breach Science 
Publishers, 1978. 

[5] E.B. Christoffel, Tiber die Gaussische Quadratur und eine Verallgemeinerung derselben, 

J. Reine Angew. Math. 55 (1858), 61-82. 
[6] J. J. Duistermaat and F. A. Griinbaum, Differential equations in the spectral parameter. 

Comm. Math. Phys. 103 (1986), 177-240. 
[7] A.J. Duran, The Stieltjes moments problem for rapidly decreasing functions, Proc. Amer. 

Math. Soc. 107 (1989), 731-741. 
[8] A.J. Duran, Orthogonal polynomials satisfying higher order difference equations, to ap- 
pear in Constr. Approx. 
[9] A.J. Duran, Symmetries for Casorati determinants of classical discrete orthogonal poly- 
nomials, to appear in Proc. Amer. Math. Soc. 
[10] A.J. Duran, Orthogonality properties for polynomials satisfying higher order difference 

equations. In preparation. 
[11] A. Erdclyi, W. Magnus, F. Oberhettinger, F. G. Tricomi, (Bateman project). Higher 

Trascendental Functions, Volume II McGraw Hill, New York, 1953. 
[12] F.A. Griinbaum and L. Haine, Orthogonal polynomials satisfying differential equations: 
the role of the Darboux transformation, in: D. Levi, L. Vinet, P. Winternitz (Eds.), 
Symmetries an Integrability of Differential Equations, CRM Proc. Lecture Notes, vol. 9, 
Amer. Math. Soc. Providence, RI, 1996, 143-154. 
[13] F. A. Griinbaum and L. Haine, Bispectral Darboux transformations: an extension of the 

Krall polynomials, Int. Math. Res. Not. 8 (1997), 359-392. 
[14] F. A. Griinbaum, L. Haine and E. Horozov, Some functions that generalize the Krall- 

Laguerre polynomials, J. Comput. Appl. Math. 106, 271-297 (1999). 
[15] F. A. Griinbaum and M. Yakimov, Discrete bispectral Darboux transformations from 
Jacobi operators. Pac. J. Math., 204, 395-431 (2002). 



X>-OPERATORS, ORTHOGONAL POLYNOMIALS AND DIFFERENCE EQUATIONS 43 



[16] J. J. Guadalupe, M. Perez, F. J. Ruiz and J. L. Varona, Asymptotic behaviour of or- 
thogonal polynomials relative to measures with mass points, Mathematika 40 (1993), 

331-344, 

[17] L. Hainc, Beyond the classical orthogonal polynomials, in; The bispcctral problem (Mon- 
treal, PQ, 1997), 47-65, CRM Proc. Lecture Notes, 14, Amer. Math. Soc, Providence, 
RI, 1998. 

[18] P. Iliev, Krall-Ja<;obi commutative algebras of partial differential operators, J. Math. 

Pures Appl. 96 (2011) 446-461. 
[19] P. Iliev, Krall-Laguerre commutative algebras of ordinary differential operators, Ann. 

Mat. Pur. Appl. DOI 10.1007/sl0231-011-0219-9 
[20] J. Kockock and R. Kockock, On a differential equation for Kooriiwinder's generalized 

Laguerre polynomials, Proc. Amor. Math. Soc. 112 (1991), 1045-1054. 
[21] J. Koekoek and R. Koekock, Differential equations for generalized Jacobi polynomials, J. 

Compt. Appl. Math. 126 (2000), 1-31. 
[22] R. Koekoek, Differential Equations for Symmetric Generalized Ultraspherical Polynomi- 
als, Trans. Amer. Math. Soc. 345 (1994), 47-72. 
[23] R. Koekoek, P. A. Lesky and L.F. Swarttouw, Hypergeometric orthogonal polynomials 

and their g-analogues, Springer Verlag, Berlin, 2008. 
[24] T. H. Koornwinder, Orthogonal polynomials with weight function (1 — x)°'(l + x)^ + 

MS{x + 1) + NS(x - 1), Canad. Math. Bull. 27 (1984), 205-214. 
[25] H. L. Krall, Certain differential equations for Tchcbycheff polynomials, Duke Math. 4 

(1938), 705-718. 

[26] H. L. Krall, On orthogonal polynomials satisfying a certain fourth order differential equa- 
tion. The Pennsylvania State College Studies, No. 6, 1940. 

[27] A. M. Krall and L. L. Littlejohn, On the classification of differential equations having 
orthogonal polynomial solutions. II, Ann. Mat. Pura Appl. 149 (1987), 77-102. 

[28] K.H. Kwon and D.W. Lee, Characterizations of Bochner-Krall orthogonal polynomials of 
Jacobi type, Constr. Appr. 19 (2003), 599-619. 

[29] K.H. Kwon, L.L. Littlejohn and G.J. Yoon, Orthogonal polynomial solutions of spectral 
type differential equations: Magnus' Conjecture, J. Appr. Th. 112 (2001), 189-215. 

[30] K.H. Kwon, L.L. Littlejohn and G.J. Yoon, Construction of differential operators hav- 
ing Bochner-Krall orthogonal polynomials as eigenfunctions, J. Math. Anal. Appl. 324 
(2006), 285-303. 

[31] O. E. Lancaster, Orthogonal Polynomials Defined by Difference Equations, Am. J. Math., 
63 (1941), 185-207. 

[32] L. L. Littlejohn, The Krall polynomials: a new class of orthogonal polynomials, Quaest. 
Math. 5 (1982), 255-265. 

[33] L. L. Littlejohn, On the classification of differential equations having orthogonal polyno- 
mial solutions, Ann. Mat. Pura Appl. 93 (1984), 35-53. 

[34] L. L. Littlejohn, An application of a new theorem on orthogonal polynomials and differ- 
ential equations, Quaest. Math. 10 (1986), 49-61. 

[35] A.F. Nikiforov, S.K. Suslov and V.B. Uvarov, Classical orthogonal polynomials of a dis- 
crete variable. Springer Verlag, Berlin, 1991. 

[36] G. Szego, Orthogonal Polynomials, American Mathematical Society, Providence, RI, 
1959. 

[37] G.J. Yoon, Darboux Transforms and orthogonal polynomials. Bull. Korean Math. Soc. 
39, (2002) 359-376. 

[38] A. Zhedanov, A method of constructing Krall's polynomials, J. Compt. Appl. Math. 107 
(1999), 1-20. 

A. J. DuRAN, Departamento de Analisis Matematico, Universidad de Sevilla, Apdo (p. 
O. BOX) 1160, 41080 Sevilla. Spain. 
E-mail address: duranSus.es 



